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Abstract 



B®-B® oscillations are observed in "self-tagged" samples of partially recon- 
structed B mesons decaying into a lepton and a charmed meson collected 
in pp collisions at ^/s = 1.8 TeV. A flavor tagging technique is employed 
which relies upon the correlation between the flavor of B mesons and the 
charge of nearby particles. We measure the flavor oscillation frequency to be 
Arn-d = 0.471^0068 ^ 0.034 ps _1 . The tagging method is also demonstrated 
in exclusive samples of — > J/ipK + and B^ — ► J /ipK*° (892), where sim- 
ilar flavor-charge correlations are observed. The tagging characteristics of 
the various samples are compared with each other, and with Monte Carlo 
simulations. 

PACS numbers: 14.40.Nd, 13.20He, 13.25Hw 
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I. INTRODUCTION 



The study of B mesons has been important for understanding the relationships between 
the weak interaction and the mass eigenstates of quarks, described in the Standard Model by 
the Cabibbo-Kobayashi-Maskawa (CKM) matrix [Q. Early studies were based on branching 
fraction and lifetime measurements. However, since the observations of B°-B° mixing, first 
in an unresolved mixture of B® and B° s by UAl ||, and then specifically for the B^ by 
ARGUS 0, a new window on the CKM matrix was opened. B° mixing, analogous to K° 
mixing, is possible via higher order weak interactions, and is governed by the mass difference 
Am between the two mass eigenstates. Unlike the K° system, the B° mixing amplitude is 
dominated by the exchange of virtual top quarks, and so provides a view of weak charged 
current transitions between a top quark and the quarks composing the B® s . 

Mixing measurements are predicated upon identifying the "flavor" of the B° meson at 
its time of formation and again when it decays, where by "flavor" we mean whether the 
meson contained a b or b quark. Determination of the initial flavor is the primary difficulty, 
as knowledge of the decay flavor is usually a byproduct of the B reconstruction, even if it is 
only partial. 

The effective size of flavor tagged B samples is a critical limitation of current measure- 
ments, especially for exclusive B reconstructions. This fact has motivated efforts to develop 
a variety of tagging techniques to fully exploit existing data. There has been considerable 
progress in recent years in utilizing a variety of tagging methods and B° samples, as illus- 
trated by the diversity of mixing measurements [[| . Even though a new generation of high 
statistics B experiments will soon come on-line ||, many tagging-based studies — such as 
CP violation in B mesons — will still be statistics limited. Thus, improvements in tagging 
capabilities will be valuable in the next generation of B experiments as well as for the current 
ones. 

We have reported in an earlier Letter || the development and application of a "self- 
tagging" method based on the proposal J7J that the electric charge of particles produced 
"near" the reconstructed B meson can be used to determine its initial flavor. Such corre- 
lations, first observed in e + e~ — > Z° — > bb events by OPAL ||, are expected to arise from 
particles produced from decays of the orbitally excited B** mesons, as well as from the 
fragmentation chain that formed the B. We refer to this approach as "Same Side Tagging" 
(SST), in contrast to other tagging methods which rely upon the other 6-hadron in the event. 

We applied SST to a large sample of B u ^ — > ID^X decays: the expected time depen- 
dent flavor oscillation was observed, and its frequency Am^ was measured with a precision 
similar to other single tagging results. In addition to the intrinsic interest of obtaining a 
supplementary measurement of Am d , this result also demonstrated that this type of tagging 
method is effective even in the complex environment of a hadron collider. A variant of this 
approach has also been studied by ALEPH |Sj in exclusively reconstructed B's at the Z° 
pole. 

In this paper we describe in detail the SST method we have developed and its previously 
reported application to B u d — > £D^*'X decays. Experimental complications surrounding 
the use of these decays are described in detail, i.e., both the cross-talk between B+ and B%, 
and the contamination from tagging on B decay products. The value of Am^, as well as the 
purity of the flavor-charge correlations, are reported. 
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This paper extends the application of SST to two fully reconstructed B decays which 
offer another test of its effectiveness: B+ -> J/ipK + and B° d -> J/^^*°(892).[| Although 
our samples are too small to yield precise tagging results, they are the largest currently in 
existence and serve as a prototype for tagging B® — > J/ipKg [|T0| , |TT]| , the centerpiece of future 
CP violation studies with B mesons 0. The tagging results from the J/tpK samples are 
compared to those from £D^*\ and also to Monte Carlo simulations. The simulation offers 
further insights into the behavior of this SST method. 

This paper is structured as follows. We review the relevant aspects of our detector and 
data collection in Sec. |TJ. Section |TTT| summarizes B°-B° mixing, and is followed by some 
remarks on tagging and a description of our specific SST method in Sec. [TV]. Same Side 
Tagging is applied to the £D^*> sample in Sec. 0, which includes discussion of B recon- 
struction, sample composition, proper decay time measurement and corrections, the tagging 
asymmetries, and finally extraction of Aro^ and tagging dilutions. This completes our main 
result. 

Having established the technique in £D^*\ we extend SST to the exclusive J/ipK modes 
in Sec. |VJ. We discuss the sample selection, the fitting method, and the resultant tagging 
dilutions. Special attention is given to handling tagging biases. Finally in Sec. |VTJ| we 



present some checks of our measurements and compare the behavior of this tagger in these 
two different types of B decays. Aspects of the data are also compared to Monte Carlo 
simulations, and the behavior of this SST method is discussed. We close with a few remarks 
concerning future applications of this type of SST method. 

II. THE CDF DETECTOR AND DATA COLLECTION 

A. Apparatus 

The data discussed here were collected using the CDF detector in the Tevatron Run 
I period during 1992-1996, and comprise approximately 110 pb _1 of integrated luminosity 
of pp collisions at yfs = 1.8 TeV. Details of the CDF detector have been previously pub- 
lished pT2"| , |l3|| , and only the features relevant to this analysis are reviewed here: the tracking 
system by which charged particles are identified and their momenta precisely measured, the 
central calorimeters for electron identification, and the muon chambers for muon identifi- 
cation. Our coordinate system is such that the (spherical) polar angle 9 is measured from 
the outgoing proton direction (+,2-axis) and the azimuthal angle </> from the plane of the 
Tevatron. 

The tracking system consists of three detectors immersed in a 1.4 T magnetic field 
generated by a superconducting solenoid 1.5 m in radius. The innermost tracking device 
is a silicon microstrip vertex detector (SVX) fl3f , which provides spatial measurements 
projected onto the plane transverse to the beam line. The SVX active region is 51 cm long 
and composed of two 25 cm long cylindrical barrels. Each barrel has four layers of silicon 



1 Reference to a specific particle state implies the charge conjugate state as well; exceptions are 
clear from the context. 
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strip detectors, ranging in radius r from 3.0 to 7.9 cm from the beam line. The impact 
parameter resolution of the SVX is aaipr) = (13 + 40/pr) /im, where px is the transverse 
momentum of the track relative to the beam line in GeV/c. The geometrical acceptance 
of the SVX is about 60% for the data presented here due to the ~ 30 cm RMS spread 
of the pp interactions along the beam line. Outside the SVX is a set of time projection 
chambers (VTX) which measure the position of the primary interaction vertex along the z- 
axis, and is in turn surrounded by the central tracking drift chamber (CTC). This 3 m long 
chamber radially spans the range from 0.3 to 1.3 m, and covers the pseudorapidity interval 
\r]\ < 1.1 (77 = — ln[tan(0/2)]) relative to the nominal interaction point. The 84 radial wire 
layers of the CTC are organized into nine "superlayers." Five "axial" superlayers consist 
of wires strung parallel to the beamline. Interspersed between these five are four "stereo" 
superlayers in which the wires are turned 3°; the two types of superlayers used together yield 
three-dimensional charged track reconstruction. Within each superlayer the wires are further 
organized into "cells" which are rotated 45° relative to the radial direction. This rotation 
assists the resolution of left-right ambiguities in track reconstruction. The CTC and SVX 
combined provide a transverse momentum resolution of <j Pt /pt ~ a/ (0.9pr) 2 + (6.6) 2 x 10~ 3 , 
with pt in GeV/c. 

Outside the magnet coil, and covering the pseudorapidity range of the SVX-CTC system, 
are electromagnetic (CEM) and, behind them, hadronic (CHA) calorimeters. They have a 
projective tower geometry with a segmentation of A</> x Arj = 15° x 0.11. The CEM is a 
lead-scintillator stack 18 radiation lengths thick. It has a resolution of 13.5%/a/Et plus a 
constant 2% added in quadrature, where E T = Esm(9), E is the measured energy of the 
cell in GeV, and 9 is its polar angle. A layer of proportional chambers (CES), embedded 
near shower maximum in the CEM, provides a more precise measurement of electromagnetic 
shower profiles both in azimuth (0) and along the beam (z) direction. The CHA is an iron- 
scintillator calorimeter 4.5 interaction lengths thick, and has a resolution of 50%/v^Et plus 
a constant 3% added in quadrature. 

The calorimeters also act as a hadron absorber for the muon chambers which surround 
them. The central muon system (CMU), consisting of four layers of drift chambers covering 
\r)\ < 0.6, can be reached by muons with px in excess of ~1.4 GeV/c. These are followed by 
60 cm of additional steel and another four layers of chambers referred to as the central muon 
upgrade (CMP). The central muon extension (CMX) covers approximately 71% of the solid 
angle for 0.6 < \r)\ < 1.0 with four free-standing conical arches composed of drift chambers 
sandwiched between scintillator (for triggering). 

The data samples of interest in this paper, inclusive electrons and muons, and dimuons 
in the mass region around the J/ip, were collected using CDF's three-level trigger system. 
The first two levels are hardware triggers, and the third level is a software trigger based 
on offline reconstruction code optimized for computational speed. Different elements of the 
trigger have varying efficiency turn-on characteristics, generally dependent upon track p-r's 
or calorimeter E^s. The behavior of the trigger has been extensively studied. Since the 
analyses presented here are largely insensitive to trigger behavior, we refer the interested 
reader to Ref. [p!4| , |l5| for detailed discussion of the triggers and their performance. 
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B. Inclusive lepton data set 

The inclusive lepton data set is composed of electron and muon triggers. Electron identi- 
fication is based on energy clusters in the CEM with an associated CTC track. The principal 
single electron trigger required a Level-2 trigger threshold of 8 GeV, and an associated 
track with > 7.5 GeV/c. The offline reconstruction requires tighter matching between 
the position of the CES cluster and the associated track (i.e., r\A<p\ < 3.0 cm and \Az\ < 5.0 
cm). The CEM cluster is also required to have a shower profile consistent with an electron 
shower, i.e., a longitudinal profile with less than 4% leakage in the hadron calorimeter, and 
a lateral profile in the CEM and CES consistent with electron test beam data. 

Muon identification is based on matching CTC tracks with track segments in the muon 
chambers. The inclusive sample is based on a Level-2 trigger with a nominal pt threshold of 
7.5 GeV/c. Each muon chamber track is required to match its associated CTC track. Track 
segments in both CMU and CMP are required to reduce backgrounds. 

The inclusive lepton triggers are the dominant contribution to our sample. However, the 
offline selection does not explicitly require that these triggers be satisfied. All events with 
a lepton track of pr > 6.0 GeV/c, and passing the above identification quality cuts, may 
enter this sample. The contribution from other triggers is small, and the bulk of events with 
lepton pt below the nominal 7.5 GeV/c threshold arise when the lepton p? reconstructed 
offline is lower than that estimated by the trigger system. Finally, only lepton candidates 
using SVX tracking information are considered, so as to be able to do precision vertexing. 



C. J ftp data set 

The J/ip sample is based on a dimuon trigger. The trigger and selection on each muon 
are similar to that for the inclusive muons described above, except for a lower nominal px 



threshold of ~2 GeV/c [fig] . The CMU-CMP requirement is also relaxed: the muon candi- 
dates may be in any of the muon chambers (CMU, CMP, or CMX), and in any combination. 
The Level-3 trigger requires the presence of two oppositely charged muon candidates with 
combined invariant mass between 2.8 and 3.4 GeV/c 2 . In offline reconstruction we further 
impose tighter track matching and require px > 1-5 GeV/c for each muon. We also require 
a minimum energy deposition of 0.5 GeV for each muon in the hadron calorimeter, as ex- 
pected for a minimum ionizing particle. Again, the dimuon sample is not explicitly required 
to have passed the dimuon trigger. 

At this stage, no SVX tracking requirement is imposed, and there are about 400, 000 
J/^'s reconstructed, with a signal-to-noise of about 10:1. Only about half of these are fully 
contained within the SVX. 



III. B°-B° MIXING 

The phenomenon of B°-B° mixing, analogous to K°-K° mixing, occurs via higher order 
weak interactions. Starting with an initially pure sample of B 0l s at proper time t = 0, the 
numbers of B° and B° mesons decaying in the interval from t to t + dt are dN(t)B0^B° and 
dN{t) 

bo^bo respectively; and they are given by 
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dNms^^Np^ +cosAmt) 

at 2tq 

dt 2r K J> K J 

where tq is the average lifetime of the two neutral B meson eigenstates, and Am is the mass 
difference between them. 

To observe mixing one must experimentally determine the flavor of the neutral B me- 
son at the times of formation and decay, a process referred to as "flavor tagging." The 
flavor at decay is usually well known from the observed decay products. The initial flavor 
determination is more difficult and is discussed in the next section. 

In an experiment with no background and perfect flavor tagging and lifetime reconstruc- 
tion, the mixing frequency Am^ can be determined from the asymmetry 

Ao(t) = -f — — f — — = cos Amt. (3) 

Tt N(t) B o^ B o + Tt N(t) B o^ 

If the flavor tag correctly identifies the B° flavor at production with only a probability 
Vo, then the amplitude of the measured asymmetry A^ eas \t) is reduced by a factor V = 
2V — 1, called the "dilution," i.e., 

A parallel series of expressions may be written when tagging -B^'s, but there is no time 
dependence, so 

Af eas \t) = {2V + -l) = V + . (5) 

Tagging charged -B's can be used to infer the flavor of the other b hadron in the event, but 
in this paper it is principally of interest as a test of the tagging method. The charged and 
neutral dilutions need not be equal, and V + can not in general be used as a measure for V . 

The uncertainty on a measurement of the asymmetry A from a sample of iV (background- 
free) events is 

a\ = (1 - A 2 V 2 )/NeV 2 ~ 1/NeV 2 , (6) 

where e is the efficiency to obtain a flavor tag for the method being employed. The figure 
of merit, eV 2 , is called the "effective tagging efficiency" of the method. 



IV. FLAVOR TAGGING 
A. Tagging methods 

There is now a considerable inventory of B® mixing measurements available [|J. Most 
rely on determining the flavor of the second 6-hadron in the event to infer the initial flavor 
of the originally reconstructed B meson. Examples include lepton tagging and jet-charge 
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tagging [|R|. We refer to these as "Opposite Side Tagging" (OST) methods. Reliance on 



the opposite-side 6-hadron can have several disadvantages. 

At the Tevatron, once one B meson is produced in the central rapidity region covered by 
CDF, the second 6-flavored hadron is present only ~ 40% of the time in this region. In the 
other ~ 60% of events the second 6-hadron is unavailable for tagging. For lepton tagging, 
there is the additional inefficiency arising from the semileptonic branching ratio of the B, as 
well as the confusion from daughter charmed particles decaying to leptons. For jet-charge 
tagging, the purity of the flavor-tag decision is reduced by the presence of charged tracks 
from the proton-antiproton remnants and possible confusion with gluon (or light quark) jets. 
Finally, tagging based on OST suffers from the inevitable degradation arising from mixing of 
the second 6-flavored hadron when it is a B°. In spite of these complications, OST methods 



have proven to be powerful tagging methods in previous mixing measurements [|T^,|T^ . 

A contrasting approach is "Same Side Tagging" (SST), which ignores the second b- 
flavored hadron and instead considers flavor-charge correlations of charged particles pro- 
duced along with the B meson of interest .0 Such correlations are expected to arise from 
particles produced in the fragmentation chain and from decays of B** mesons. 

A simplified picture of the possible fragmentation paths for a b quark is displayed in 
Fig. [IJ. If the b quark combines with a u quark to form a B+, then the remaining u quark 
may combine with a d quark to form a 7r~. Alternatively, if the b quark fragments to form a 
B®, the correlated pion would be a 7r + . These correlations are the same as those produced 
from B** decays, such as B*f° — > £?i* )+ 7r _ or B** + — > I^* )0 7r + . We do not attempt to 
differentiate the sources of correlated pions. 

In this simple picture of B-n correlations, naive isospin considerations imply that the 
tagging dilutions for _B°'s and £?+ 's should be the same. However, this need not be the 



case 



|21~H , and we make no such assumption. Furthermore, we generically refer to the tagging 



particle as a pion, although we do not attempt to experimentally identify it as one. 



B. Same Side Tagging algorithm 

General considerations of correlations between B flavor and particles produced in frag- 
mentation offer only qualitative guidance in constructing an SST algorithm. String frag- 
mentation models indicate that the velocity of fragmentation particles are close to that of 
the B, and similarly for pions from B** decays. Motivated by this observation, a number 
of variables were studied for selecting a tagging track using data and Monte Carlo simula- 
tions, among them: (i) the maximum px track, (ii) the minimum S-track mass (using the 
pion mass), (Hi) the minimum AR = a/ (Ar/) 2 + (A0) 2 between the B and track, (iv) the 
minimum of the track momentum component transverse to the combined momentum of the 



Jet-charge tagging has been extended by combining the opposite and same side jet-charges in 
Z° — > bb [flSfl . A same side jet-charge tag is clearly correlated with the SST approach of this paper, 
but the philosophy is different. The jet-charge method is based on a weighted average of charged 
tracks reflecting the primary quark's charge |p0| , while the proposal of Ref. is based more on 
selecting a specific charged particle to determine the flavor. 
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B (Pb) plus track (Ptr) momentum (p™ 1 , see Fig. ||), and (v) the maximum of the track 
momentum component along the 5-track system momentum (p r £ l ), as well as several others. 
We found that these five variables have similar performance, and moreover were highly cor- 
related in selecting the same track as the tag. Future studies with higher statistics samples 
may enable one to optimize the choice, but we were unable to identify one method as clearly 
superior. We chose to use p™ 1 , as this variable was among the best for correctly identifying 
the flavor (i.e., had a large T>), and it seemed less vulnerable to tagging on B decay products 
missed in partial B reconstructions (Sec. |V hp ,p| 

For our specific SST algorithm, we consider all charged particles that pass through all 
stereo layers of the CTC and are within the r]-<f> cone of radius fractionAi? = 0.7 centered 
along the direction of the B meson. If the B is partially reconstructed, we approximate 
this direction with the momentum sum from the partial reconstruction. Tracks are required 
to be consistent with having originated from the fragmentation chain or the decay of B** 
mesons, i.e., coming from the primary vertex of the pp-interaction. This translates into the 
demand that tracks must have at least 3 out of 4 SVX hits, c?o/Vo < 3 where do is the 
distance of closest approach of the track trajectory to the primary vertex when projected 
onto the plane transverse to the beam line (r-<p plane) and ctq is the estimated error on do, 
and the closest approach in z must be within 5 cm of the primary vertex. 

Due to chamber design, the CTC is known to have a lower reconstruction efficiency for 
negative tracks compared to positive ones at low p T (Sec. [VI C 3|) . To suppress this bias, all 
candidate tracks must have a px above a threshold of pr(SST) = 400 MeV/c. 

At this point, more than one candidate tag may be available for a given B. To select the 
tag, we choose the candidate track with the smallest p™ 1 '. 

A B is tagged if there is at least one track that satisfies these selection requirements. 
The fraction of B candidates with a tag is the tagging efficiency, and it is ~ 60-70% for this 
algorithm in our data. 

V. FLAVOR OSCILLATIONS IN THE LEPTON+CHARM SAMPLE 

We apply our SST method to a sample of B u>d decays to a lepton plus charmed meson. 
We form the asymmetry, analogous to Eq. (|3p, between the decay flavor and the charge 
of the tag track, and we fit this asymmetry using a \ 2 minimization to obtain Am^. As 
a by-product, the tagging dilutions are also determined. As we are henceforth concerned 
specifically with £>+ and B%, the subscripts are suppressed for the remainder of this paper. 

The incomplete reconstruction of the B's introduces several complications: (i) missing 
decay products means that the precise /?7-factor to compute the proper decay time is not 
known; (ii) a missed charged decay product results in a B + being classified as a B° and 
vice versa; and (in) a missed charged decay product may be chosen as the tag, biasing 
the asymmetry. The latter two issues are the principal subtleties of this analysis, and 
necessitate careful consideration of the composition of the sample. Not all the branching 



3 If not for this issue, our studies tended to favor p r £ l , essentially the same variable employed in 
Ref. H for tagging exclusively reconstructed B samples. 
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ratio information required is well known, and when not, we rely internally on our data set. 
Because the unknown sample composition parameters depend themselves on other sample 
composition parameters we use an enlarged \ 2 function to fit globally for Am^ and the 
unknown composition parameters. 

We first describe the sample selection and then discuss issues of sample composition. 
The proper time measurement, and the corrections for missing particles, is fairly standard, 
but B° <-» B + cross-talk introduces additional corrections. We then discuss the measured 
and expected flavor-charge asymmetry given the complications of the sample composition, 
including the biases of tagging on B decay products. We finally discuss the \ 2 fit; results, 
and the effects of systematic uncertainties on Am^ and the tagging dilutions. 



A. B candidate selection 

We use partially reconstructed i?'s consisting of a lepton and a charmed meson. A 
particular B reconstruction does not necessarily arise from a unique sequence of bottom and 
charm decay modes when there are unidentified decay products (Sec. |V Bp . We therefore 
refer to the various B reconstructions as "decay signatures," and use the predominant decay 
sequence as a label. The samples of B 0, s consist of four decay signatures, one B° — > vl + D~ 
signature and three B° — > v£ + D*~: 

(7) 
(8) 
(9) 

(10) 

where we adopt the convention that a 7r from a D* or D** decay is labeled by a "*" or 
subscript. For the I? +, s, we use only one decay signature: 

B + -> u£+D°, D° -> K+tj-. (11) 

As noted above, the decay signatures do not represent a specific sequence of decays; they in 



B°- 


-> vt 


D-, 


D~ - 


■* K + 7l 7T 






B°- 


-> vt 


D*- 


D*- 




7J°- 


-> K + TT- 


B°- 


-> vt 


D*- 


D*~ 




7J°- 


-> K + TT- 


B°- 


-> vt 


D*- 


D*- 


-»■ D TT^ , 


7J°- 


-> K + it- 



fact include several sequences, for instance, Eq. (p~T| ) includes the decay chain B + — ► u£ + D*° 
followed by D*° — > I^ ^ and D° — > K + tt~, where the 7r° is not identified. 

The B selection starts with the inclusive lepton (e and fi) samples of Sec. |11B| . The tracks 
of the daughters must lie within a cone of AR =1.0 around the lepton, pass through 
all nine CTC superlayers, have enough hits for good track reconstruction, and satisfy a Pt 
requirement (see Table |). All tracks (except one in the case of D° — > K + 7t~7t~tt + ) must use 
SVX information, and they must also be consistent with originating in the vicinity of the 
same primary vertex. The candidate tracks must form an invariant mass in a loose window 
around the nominal D mass, where all permutations of mass assignments consistent with 
the charm hypothesis are attempted. 

The candidate tracks are combined in a fit constraining them to a D decay vertex; x 2 an d 
mass window cuts are imposed. With the D vertex established, we select the primary vertex 
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from thoseQ reconstructed in the VTX as the one nearest in z to the D. The transverse 
coordinates of the primary vertex are obtained from a z dependent average beam position, 
as measured by the SVX over a large number of collisions recorded under identical Tevatron 
operating conditions. We require the D tracks to be displaced from this primary vertex 
(do/ (Jo cut in Table |), and the projected transverse distance L xy (D) between the D vertex 
and the primary vertex to be greater than its uncertainty <j£ xy (L xy (.D)/<7£ xy cut in Table |). 
The projected distance L xy (£)) is defined as 

xy( } = [MW\ ' ( J 

where the two vertices are given by the position vectors x pr i m and xd, and the D transverse 
momentum is Pt{D). 

We next find the B vertex. For the B° — > u£ + D*~ signature the lepton and the ir~ from 
the D*~ decay both come from the B decay point. We fit for the B vertex by intersecting 
the lepton and the 7r~ tracks, and require that the D points back to it. For the B° — > v£ + D~ 
or B + — ► u£ + D° signatures there is no additional track emerging from the B vertex. The 
D is projected back to the lepton track and their intersection determines the B vertex, as 
sketched in Fig. [| A loose cut is applied to the D proper decay length relative to the B 
vertex (c£d in Table [[]). The charges of the lepton and the charm candidates are required to 
be consistent with the decay of a single B, i.e., a £ ± K ± correlation. 

The B° — > v£ + D*~ decays followed by D*~ — ► D°n~ also contribute to the £ + D° samples. 
The separation between B° and B + is improved by removing all £ + D° candidates that also 
participate in the £ + D*~ reconstruction. We define a D*~ candidate as a valid D° candidate 
with a 7t~ candidate that makes the mass difference m(D°7r~) — m(D°) consistent with 
the known mass difference between the D*~ and D° |22|. Since the m(D°7T~) — m(D°) 
distribution for real .D*~'s is very narrow (~1 MeV), this removal is very efficient once the 
tt~ is reconstructed. 

The numbers of B candidates are extracted from a fit of the charm mass distributions. 
Figure |] shows the invariant mass distributions (solid histogram) for the four channels of 
exclusively reconstructed charm. The signal components of the D mass distributions are 
modeled by Gaussians, and the combinatorial backgrounds by linear functions (solid curves 
in Fig. |). 

The dashed histograms in Fig. || represent the U D" mass distributions for B candidates 
where the lepton and the kaon have the wrong charge correlation (£ ± K =F ). These "wrong- 
sign" events can be combinatorial background, as well as cases where there was a real D and 
a fake lepton. The absence of a peak in the wrong-sign U D" mass distribution demonstrates 
that the right-sign sample is a clean signal of £D ( -* S) pairs coming from single B mesons. 

In the case of the decay £ + D*~ , D° — ► K + n~n°, the tt° is lost, and the Kir invariant 
mass distribution has a broad excess below the D mass. However, in the m(iT7r7r*) —m(Kir) 
distribution a relatively narrow peak appears at the value m(D*~) — m(D°), as seen in 
Fig. j|. We parameterize the combinatorial background by the shape of the wrong-sign 



4 It is not uncommon to have multiple pp interactions in a single bunch crossing at the higher 
Tevatron luminosities. 
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{£ K^) distribution (lower dashed curve). This shape, combined with the signal function, 
is then fit to the right-sign {£ K ) data, and is shown by the solid curve in Fig. |5[ 

This completes our sample selection, which has yielded almost 10,000 B mesons. How- 
ever, before we can use them, several other issues must be addressed. 

B. The composition of the ID^*> sample 
1. B° <-> B + cross-talk 

As noted earlier, the SST correlation depends on whether the B meson was charged 
or neutral. However, only the ground state charm mesons and one D* decay mode were 
reconstructed in the previous section, and the existence of the intermediate resonances D* 
and D** introduce contamination from B + decays into B° decay signatures, and vice versa, 
when charged decay daughters are unidentified or unreconstructed. We disentangle this 
cross-talk by relating the charged and neutral B fractions to the number of reconstructed 
charm mesons through relative branching ratios and reconstruction efficiencies. This section 
details this connection. 

There are two causes of the B° «-> B + cross-talk in this analysis: 

• Missing the n~ from the D*~ decay. For example, 

B + -> u£+D° (13) 

can be mimicked by the decay sequence 

B° -> u£ + D*~ 
D*- - D\; (14) 

if the 7r~ is not part of the reconstruction. 

• B decays to D**-mesons. The decay chain 

B° -> u£ + D**- 
D**-^D°7r- (15) 

will also mimic the £ + D° signature of the B + when the 7r~ is unidentified. 

The first case is of concern as it is not unusual for the momentum of tc~ to fall below our 
Pt cut. The 7r~ tends to be soft because of the small energy release in the decay, coupled 
with the modest boost of most of our B's. We identify the 7r~ only with some efficiency 
e(vr*). 

In the other case, we do not attempt to find the 7r#* from the D**. There are four 
expected orbitally excited D** resonances (see Table |II|), some of which decay into Dti, 
others to D*tt, and one to both. The total decay rate to these states is not well known, and 
the proportions of the four possible D** states are almost totally unknown. There is evidence 
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that the £)i(2420) and -D2(2460) states are produced at some level in B decays |[23|| . There 
may also be non-resonant D^ir production (B — > isiD^n) [23], which has the same cross- 
talk effect. It would be extraordinarily difficult to distinguish these decays from the two 
D** resonances which are predicted to be broad by Heavy Quark Effective Theory [p4j| . We 
therefore subsume the effects of four D** resonances, as well as the four-body semileptonic 
decay of the B meson, into our treatment of "_D**" 's. 

The complete picture of the decay chains is more complicated, since both B° and B + 
decay into "D**"'s, and D**~ and D**° decay into both D*~ and D*°, as well as D~ and 
D°. The full complexity of the decays is illustrated in the state diagram shown in Fig. ^|. To 
reiterate our terminology, a specific sequence of decays in Fig. ^| is called a "decay chain, " and 
the reconstructed final state is a "decay signature. " Several decay chains may contribute 
to the same decay signature. Decay chains in which the B decays directly into a decay 
signature (i.e., no particles except the neutrino are missed) are called "direct decay chains." 
Equations (p"3"D, (PI), and (p~5| ) are examples of decay chains; Eq. flUf ) is also a direct decay 
chain. Each of the five decay signatures considered in this analysis consists of several decay 
chains: three for every £ + D*~ , nine for the £ + D~ , and twelve for the £ + D°. 



2. Determining the sample composition 

Due to the B° <-> B + cross-talk, a simple computation of the time-dependent charge- 
flavor asymmetry of Eq. (||) for a £D^ decay signature will result in a weighted average of the 
B° and B + asymmetries [Eqs. (^) and (^)j. The weighting is determined by the fractional 
contributions of B° and B + decays to that decay signature; we call these fractions the 
"sample composition." The fraction of B° and B + decays in a decay signature is essentially 
determined by the branching ratios and reconstruction efficiencies for each decay chain 
contributing to that signature. Since only fractions are involved, it is convenient to use 
ratios of branching ratios and relative efficiencies. These quantities, along with the B° and 
B + lifetimes, fully describe the sample composition as a function of proper time and are 
referred to as the "sample composition parameters." We now discuss how we determine 
the fractions of B° and B + mesons contributing to a signature, given our choice of sample 
composition parameters. 

We tabulate all possible decay chains that feed into each signature, and classify them as 
originating from a B° or B + . For compactness, we label decay signatures by k, and decay 
chains by h. The symbol U B° —>■ h G k" is interpreted as "the decay chain h originates 
from a B° and contributes to the decay signature k." We write the fraction of neutral and 
charged mesons contributing to a decay signature k as 



WW 

i{N° k (t)+K(t)Y 



Kit) ee . r * */' ., , (17) 



where the dN°' + (t) are the numbers of events of signature k originating from B° or B + 
which decayed in the proper time interval from t to t + dt. These numbers are sums over all 
the decay chains h resulting in the signature k: 
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^ = ^ £ *<* - - A> (18) 

«m = » s g!i E e(B ^ ft)£(B + _ h) , (19) 

+ B+^h&k 

where we assume equal numbers of B°'s and 5 +, s are produced (i.e., Nb = 2N + = 2N°); 
r 0i+ are the B 0,+ lifetimes, and B(B 0,+ —>■ fa G fc) and e(i?°' + — > fa 6 fc) are the branching 
ratios and reconstruction efficiencies^ of a decaying through the chain fa and resulting 
in signature k. The sums for the two mesons are different since they are over a different set 
of decay chains for signature k. Knowing all the branching ratios and efficiencies, we can 
calculate the sample fractions '. 

The efficiencies e(B 0,+ — > fa G k) share common factors across decay chains. Since only 
the ratios are needed in Eqs. ([16]) and fllTD , we express the efficiencies relative to the direct 
decay chain d, for signature k, 



D _ e(£W ^hek) 
s kh = e ( 5 W ^dek)' 



(20) 



The superscripts "(fa)" and "(d)" represent the charge of the B which originated the fa and 
d chains; and the U D" superscript is a reminder that these relative efficiencies are largely 
determined by the type of charm mesons in the decay chain. For the direct decay chain 
fa = d, and so e® d = 1. We determine e® h for each decay chain fa from a Monte Carlo 
simulation as discussed in the next section. 

Similarly, only the branching ratios relative to the the direct chain branching ratio are 
required here, i.e., B(B^ — > fa G k)/B(B^ d ' — > d G k). Rather than attempting to use each 
branching ratio explicitly, not all of which are well known, we can re-express the required ra- 
tios in terms of a few relative branching fractions by using a few simplifying assumptions. We 
outline this process by considering a specific example using the B° signature k = U £ + D*~. 17 
The direct chain is B° — > i/£ + D*~, and there are two "/)**" chains ( 7r** s are unidentified) : 
B° -> v£ + D**- followed by D**~ -> D*'^, and B + -> is£+D**° with D**° £>*~7r+. If 
we index these three chains by d, a and b respectively, then the branching ratios relative to 
the direct chain are: 

B{B° -4 fl £t) _ B(B° -> v£+D**-)B(D**- -> g^jQ 

BjB^TJ^k) ~ £(5° -»• v£+D*-) ^ ' 

and 

£(£+ — » b G fc) _ g(g+ -> v£ + D**°)B(D**° -> D*-7T+) 

£(5° ^dek) ~ £(5° -»• u£+D*-) ' ^ ' 

The ratio of semileptonic I? decays can conveniently be re-expressed using ratios rela- 
tive to the inclusive branching fraction to the lowest-lying D state, including decays via 
intermediate D* and D** states, B(B -> v£DX): 



5 We apply the term "reconstruction" only to those parts of the decay we identify; neutrinos and 
decay products which are not part of the decay signature are not included. 
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J0 ' + B(B°>+^v£DX) { } 

_ B(B 0>+ -> vlD*) 
J o,+ B(B°>+^v£DX) y 1 

= B(B°' + -> uiD**) 
/o ' + ~ B(B°'+ -> ^DI) ' ( j 

We assume that all the charged and neutral ratios are equal, e.g., f = f = /+. Since 
the D* and D** decay strongly they all ultimately result in a DX signature, and thus 
f + f* + /** = 1. Because the B — > v£D** fractions are the least well known, we elect as 
our two independent parameters: 

Rf = f*/f = 2.5 ±0.6, (26) 
/** = l-f-f* = 0.36 ±0.12, (27) 



where the values are derived from world averages [^2J and the f+f* average from CLEO |25 
We may now express Eq. (|21~D as 



B{B° ^aek) f** B(D **- D *- o ) (28) 



solely in terms of /'s and D branching ratios. On the other hand, in Eq. (P2|) we can use 
B(B + -> u£ + D**°) /** ,6(5+ -> i^+I3 X) 



£(5° -> u£+D*-) f* B(B° -> u£+D-X) ' 



(29) 



where the ratio of the inclusive branching fractions to semileptonic decays of _B + relative to 
.B must be taken into account. 

The ratio B(B + — > vP. + D°X) / B{B Q — > z/£ + .D~X) can be approximated by the ratio of 
the _B + and _B° inclusive semileptonic branching ratios B s i(B + ) / B s i(B°) . According to the 
spectator model, the semileptonic width T s i is expected to be the same for the B° and B + . 
The ratio of the semileptonic branching ratios (B s i = T s i/T tot ) for the B + and B° is then 
proportional to the ratio of their lifetimes, i.e., 

Bsi(B+) = T sl (B+)/T tot (B+) = Ttot(B°) = t+ 

B sl {B») T sl {B»)/T tot {B») T tot {B+) r ' 1 ' 

This allows us to also rewrite Eq. (f22|) in terms of /'s and D branching ratios as 



s^7^) = ^7 7 ( J ' (31) 

with the I? lifetimes as the only additional parameter. We use the 5° lifetime and the ratio 
of B lifetimes as our two independent parameters, with the values 

cr = 468 ± 18 /mi (32) 
= 1.02 ±0.05 (33) 

TO 
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obtained from world averages p|] . 



The final branching ratios required are those for the charm decays. For B(D** — > P*7r«), 
we need the fraction of D** states that decay via P*7r« or Dir**. Isospin symmetry gives 
relative exclusive branching ratios for a particular D** species decaying to a D or D*, such 

as 

B(D**- -> £(*)-7T°J = B(D**° -> ggVgJ = 1 

Bp*- -»■ D(*)°ir**) Bp** -> DW-7T+) 2' 1 j 

As noted before, we use the symbol "£)**" to represent the sum over all four D** states 
(Table ^TJ) as well as two non-resonant channels. The various "/)**" states, however, decay 
differently to D and D*. Reference to a B — > P** — > P/*) decay chain implies summing over 
all possible "_D**" states. We use Py to denote the inclusive probability that a P** decay 
yields a P* as opposed to a D, and it is given by 



B(B — > P** -> P*vr„ 

v 



Bp -> P** -> P*vr„) + Bp -> P** -> Pvr«) 



(35) 



also depends upon the relative fractions of the various -B — ► D** decays since is an 
effective average over all the D** states. 
Equation (p8|) then becomes 

0>W, (36 ) 



B(B° ^dek) f* \3 

where the "1/3" is the isospin factor [similar to Eq. fl34|)1. A parallel expression may be 
written for Eq. (|3~I|). Py is poorly known and is often assumed to lie between 0.34 and 



0.78 [26|. However, it can be (weakly) constrained by our data, and we therefore let it vary 
as a free fit parameter in our Am d fit (Section [VF|). 

For the £ + D~ and £ + D° decay signatures, we also need D* branching fractions. The D*° 
always decays to a D° with a 7r° or photon, and the signature is always £ + D°. On the other 
hand the D*~ has two decay channels which feed into different signatures. These ratios are 



well known 22 : 



B{D*- -> D\;) = (68.3 ± 1.3)% (37) 
Bp*- -> p-vr°) = (31.7 ± 0.8)%. (38) 

From this £ + D*~ example we have seen the basic ingredients for determining the sample 
composition. In order to use a general notation, we define the relative ratio: 

_ BpW -> h e k) r d 
(Pkh - B{B^ ^v£DX)r h ' { ' 

where by B^ we denote the B charge state for the direct decay chain, and by its lifetime. 
Th is the lifetime of the B^ from which the chain h originates. In Eq. ( |39"|) the Td/rh ratio is 
included in order to cancel out the lifetime ratio that may appear in the branching ratios B 
by Eq. (|30D (e.g., in Eq. (|3l|) ) so that the 4>khS depend only on branching ratios averaged over 
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both B meson species. The <pkh^> are compiled in Table |Tl| along with the 7r* reconstruction 
efficiency factor, which is discussed in the next section. 

We then determine the sample composition fractions J^. ,+ {t) for signature k from the 
numbers of B°> + mesons dN°> + (t) [Eqs. (H) and (JlJD] as, 

d -^ = M kd e-^ J2 ( 4 °) 

B+^h&k 

with the normalization factor 

__ N B B(BW -> u£DX)e(B^ -> d e fc) 

•M^ = • (42) 

2^ 

When calculating the ratios for this factor cancels out. It is with the F^it) fractions 

that we fully quantify the sample composition. 



3. Reconstruction efficiencies 

We use a Monte Carlo simulation to calculate the relative reconstruction efficiencies 
e^ h for each decay chain h contributing to signature k relative to the direct decay chain 
for k. Many systematic effects cancel out in these ratios of lepton+charm reconstruction 
efficiencies. In fact, these ratios depend almost exclusively on the decay kinematics, which 
are reliably simulated. 

We use our single B Monte Carlo generator (App. [A 1| ) to produce and decay B mesons, 
and we pass them through the standard CDF detector simulation. We then apply the 
selection cuts and calculate the relative efficiencies. The e® h vary from about 0.2 to 1.5, 
with most of the variation arising from the effects of the fixed lepton pj- threshold on the 



reconstruction of the various D states p7 |. 



One last efficiency is needed. The D*~ — > D°tt~ reconstruction efficiency includes a 
contribution for the efficiency of finding the 7r~, which cancels out in the ratios e^ h . Loss of 
the 7r~ makes D*~ — > D tt~ look like a D°. Since D°7r~ candidates are removed from the 
D° sample, we need to know the absolute efficiency e(7r*) to quantify the separation of the 
D*~ and D° signatures. 

We use data rather than Monte Carlo to determine e(7T*), since the absolute detector 
response for such low pt particles is difficult to model. We use a related quantity, which 
can both be calculated from e{n*) and other sample composition parameters, as well as be 
measured in data. This quantity is the fraction of l + D*~ candidates reconstructed out of 
the entire D° — > K + tt~' sample (i.e., before the D* removal), 

N{£+D°(7r±)y y } 

where £ + D°(7vf) si gnifies £ + D° candidates before D* removal. The measurement in data, 
jj*{meas) ^ j g accc ,mplished by fitting the £ + D*~ and £ + D° (without D*~ removal) mass dis- 
tributions simultaneously, and returns R*( meas ) = 0.249 ± 0.008. 
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The calculation of R* consists of summing over all the decay chains which give the desired 
signatures. Each term is weighted by reconstruction efficiencies. The denominator sums over 
all decay chains which have a D° in their final states, including D*~ decays: 



J 



e+D Q x 



(t)]dt 



(N B e(u£ + D°)B{B + -> ut+D°X) 
{ 2^ 



+ T ^ kh€ kh 
B°-^h£e+D° 



(44) 



where the lifetime factors result from integrating the exponential factors over time. 

The numerator of Eq. (f4"5|), on the other hand, only sums over those decay chains which 
give a D*~, and is then multiplied by e(7r*) to make it the number of D*~ J s which are fully 
reconstructed, i.e., 



N(£ + D*- 



j J 

N B e(v£ + D°)B{B+ -> u£+D°X) 
2^7 



(t)]dt 



x 

x B(D 



r + (r^P y e a D ) + To(f*e° + /**^Pye?) 



(45) 



We have explicitly substituted the sample composition parameters <?W S from Table |T| 
in the square bracket term since it is relatively simple in this case. The subscripts on 
the relative efficiencies refer to the following chains: (a) B + — > £ + D**°, D**° — > D*~7r+; 
(6) P° -> and (c) B° -> £+0**-, £>**- -> D*~7r^. All D*~'s decay to D°tt-. We 

see that the ratio of these two expressions, the prediction for R*, is directly proportional to 
e(7r*), and only depends upon previously defined sample composition parameters. When the 
sample composition dependent prediction for R* is constrained to the value R*( meas ) in the 
Am d fit, we find that e(vr*) « 0.85 (Sec. |VF2D . 



^. Summary of the sample composition 

The fractions of the 5° and B + decays in each of the five decay signatures are described 
by a set of sample composition parameters. Among them, Rf, /**, and t + /tq are obtained 
from other experiments, and the e® h are calculated from Monte Carlo simulation. The 
parameter e(7r*) is expressed in terms of the other sample composition parameters (via R*) 
and R*( meas ) (obtained from the data). The final parameter, Py, will be a free parameter to 
be determined in the Am^ fit. 
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Measuring B -B oscillations also requires the determination of the proper time of the 
B decay. This will be discussed next, but sample composition effects must be included there 
as well. 



C. Proper time of the B decay 

The true proper time t of a B decay may be determined by using its projected transverse 
decay length relative to the primary vertex L xy (B) (following Eq. ([12]) ), by 

m{B) 
'MB)' 



^=M£)^, (46) 



where m(B) is the mass of the B and Pt{B) is its transverse momentum. Since the B is only 
partially reconstructed here, we use Monte Carlo-derived average corrections JC k h relating 
the reconstructed parts of the transverse momentum pt{£D) to that of the complete pt{B), 
i.e., 

^i^m (47) 



Pt 



(B)/ 



MC 



for decay chain h contributing to signature k. 

The /C-factors are determined from the same simulation (App. |A 1|) as the efficiencies e kh . 
An example of a pt{£D) /pt(B) distribution is shown in Fig. [7|. The distribution is relatively 
well concentrated because the lepton trigger threshold favors decays where the neutrino takes 
only a small portion of pt{B), thereby making the £D system a fair representation of the B. 
The direct decay chains have means of ~85%, and RMS's of ~ 10%. Also shown in Fig. [7] 
is the mean of pt{£D)/pt{B) as a function of the £D mass m(£D); less of a correction is 
needed the closer m(£D) is to the B mass. We improve our resolution by using a m(£D) 
dependent correction on an event-by-event basis. 

The correction factor varies with decay chain, so the complete scale factor, /C fc ' + , for 
signature A; is a sample composition- weighted average of the /Cfch's, 

for _B°'s, and an analogous expression for /C^ . In order to simplify averaging over the 
sample composition and cancel systematic uncertainties, we replace K, k h in Eq. ([48]) by 
K-kd x (K-kh/fc-kd), where d is the direct chain contributing to k. We factor K. kd outside the 
summation leaving the ratio K, k h/K, k( i behind. The set of factors we then use are the K, k d 
with the m(£D) dependent corrections, and the /C^/ZC^ averaged over m(£D) (where the 
m(£D) dependence largely cancels out). 

The factors JC k and JC k are different by virtue of the summation over different decay 
chains for S°'s and B +, s. The dependence of the sample composition T k on the lifetimes is 
accounted for by using the corrected times, 

ctr^L^^n^ct (49) 
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as an estimate of the true proper time t in the sample composition fractions, e.g., for Eq. ( |i~6|) 
we write 

■F°(rt° rt + ) = ft N k( ct k) /-on 

The use of pt{£D) rather than the true Pt{B) smears the ct distribution in addition to 
the average shift considered above. The difference between the reconstructed proper decay 
distance ct°^ + and the true distance ci is (suppressing most super- and subscripts) 

Act = ct°' + - d = A(L xy (5)//3 TT ). (51) 

Approximating 1//?t7 with its mean value (l/^Tl) gives 

which illustrates the effect of the reconstruction resolution via the AL xy (B) term, and the 
additional smearing due to our average corrections by A(1//?t7)- 

An example of the simulated AL xy (B) distribution is shown in Fig. || It has a Gaussian- 
like shape and an average resolution of a few hundred microns. Also shown is the fractional 
PtJ distribution, which is sharply peaked (RMS ~14%), and is essentially a mirror image 
of PjP/pt (Fig- 0)- The combined effect of both factors is shown by the Act distribution in 
Fig. |, it has an RMS of 140/im. 

Given the linear dependence of Act on the proper time in Eq. (|5*2]), we parameterize the 
ct resolution as 

a ct (ct) =a ct (0)+bx ct. (53) 

We use the RMS spread of the Act distribution for bin ct as the resolution a ct (ct), and 
fit the RMS values of the various ct bins for the slope and offset of Eq. ([53]). The linear 
model works well as seen for the sample chain shown in Fig. 13. This process is repeated 



for all five direct decay chains, and the results are listed in Table [IV]. Each chain has a 
somewhat different slope, but the intercepts are similar to the intrinsic detector resolution 
of ~ 40-50 /iin obtained when vertexing pairs of high px tracks at low ct [[15[ . 

The different decay chains that compose a decay signature are topologically similar. 
Simulation shows that the ct dependence of the ct resolution among the decay chains within 
a signature are very similar. We make the (small) sample composition correction to the ct 
resolution for signature k by approximating it as 



°f{ctkd) = P^o&(ct u ), (54) 

[Ckd/ct 

where cr^( c ^) is the parameterization of Eq. ([53]) for the direct chain d, and the bar 

indicates an average over contributing chains while the angle brackets denote an average 

over ct. Thus {(Jkd)ct is the ct-averaged ct resolution for direct chain d, and (o"fc) ct is the 

sample composition-weighted average, over all decay chains contributing to k, of the ct- 

averaged resolution. The parameter {<Jk) ct not only reflects the different ct resolutions of the 

various decay chains, but also the fact that the earlier use of the average correction factor 
-o, 



/C fc ' , rather than JC kh , introduces additional smearing |27[| . 
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D. Tagging tD^*> and the sample composition 



We apply the SST algorithm to the £D^ sample and find that ~70% of the events are 
tagged. We classify events for each decay signature as having the "unmixed" lepton-tag 
charge combination (e.g., £ + tt + for B 0, s and £ + n~ for B +, s), or the "mixed" one with the 
inverted 7r charge. Each set is further subdivided into 6 ct kd bins,[] where ct kd is the proper 
time obtained using the direct decay chain d correction factor for signature k (like Eq. ([49D, 
but only using fCkd)- 

The charm mass distribution for each of these ctkd subsamples is fit to a Gaussian signal 
plus linear background. The mean and width of the Gaussian, and the background slope, 
are all constrained to the same value for all the subsamples of a given signature. The fitted 
numbers of unmixed (Nj?(ct kd )) and mixed (Nj^ (ct kd )) events for signature k in the discrete 
ctkd bin are then used to compute the measured asymmetry, 

Ameas),. s. _ (ct kd ) ~ (ct k d) 

Ak [Ctkd) = NV(ct kd ) + Nl;i(ct kd y (55) 

Numerically, the value of the ct kd bin center is chosen as the average over the candidates' 
ct kd s in the bin, thus accounting for the nonuniform ct distribution from the exponential 
decay. 

Denoting the true asymmetries for B° and B + as A and A + , one has for a pure, perfectly 
identified B° sample the "predicted" asymmetry A k [o](ct) = A°(ct), where "fc[0]" indicates 
that k is a B° signature. When k also includes B + decay chains, one has 

A m (ct kd ) = F° k (ct° k , ctt)A°(ctl) + F+(ctt, ct° k ){-A + (ctt)}, (56) 

for the prediction. The true asymmetries are combined in a sample composition-weighted 
average, with the fractional contributions J-Z ,+ from Eqs. ([TB]) and (|17]). Furthermore, the 
observed proper time must now be corrected for the sample composition by using the ct° k ,+ 



from Eq. (f49|) . The A + term appears with a negative sign since the charge of the flavor- 
correlated tag is reversed when a B + is mistaken for a. B°. A similar expression, albeit with 
signs reversed, holds for a B + signature. 

A further correction for A k is necessary because there is the possibility of selecting the 
7r^ from a D** decay [see Eqn. (|15|)1 as a tag by mistake. No attempt was made in the sample 
selection to identify vr^'s. The lepton and a 7r^. tag almost alwaysf] have the right charge 
correlation for an unmixed B, given the apparent charge of the reconstructed B. The tt^ 
contribution is quantified by the relative number of D** y s present (/**) and the probability 
£ of selecting the 7T^, 8jS cL tag in a tagged event in which a 7r^ was produced. With this 
definition of £ the effect of the tagging algorithm is separated from the D** branching ratios. 



We cannot use ct k ' to bin the data because the sample composition is not completely known 
until after the binned data have been fit (Sec |VF ). 

7 The B + chain that cascades through D**° — > D*~ — > D°, tags on the 7r+ , and loses the 7r~ in 
the reconstruction, is an exception. However, this has a small contribution, i.e., JT+>** ~o.5% in 
Eq. fl57|). 
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We split the B° and B + decays into those with and without a 7r^, and define '** 
as the fraction of decay signature k in which a 7T^ was produced. jF°'** is calculated in 
the same way as in Eq. ([16]), except that here the numerator is a sum only of the decay 
chains involving a nf^ from B°. J^'** is calculated analogously. Only a fraction £ of 7r^'s are 
selected as tags, and we split the B° components into — and and similarly 

for B +, s. We then generalize Eq. ([56]) to include 7r^ tags in the prediction for the measured 
asymmetry: 

Am = W - e^**M° + e^'**{-i} 

+ W - + ^'"{+1}, (57) 

where the —1 (+1) asymmetry factors in the second (fourth) term reflect the perfect corre- 
lation of nf^ tags. 

All relevant effects for a mixing measurement using SST are contained in Eq. ([5j]); it 
describes the observed asymmetry Ak[o] given the true asymmetries A°' + , the nf^ tagging 
probability £, and the sample composition jF's. 

E. Determination of the tt^, fraction £ 

The 7r^ tagging probability £ depends on the tagging algorithm, the kinematics and 
geometry of the B and D** decays, as well as the characteristics of the fragmentation and 
underlying event tracks. We use a full event simulation (App. |A 2 a| ) to model the decay 
kinematics and geometry — which it does reliably — to obtain the ct dependence of £, denoted 
by £mc( c ^)- The decay kinematics and geometry determine the £ shape, whereas the relative 
competition between the nf^ and the other tracks to become the tag affects the overall 7n^ 
tagging probability. This observation enables us to use our data to determine the global 
normalization of £, instead of relying on the simulation's modeling of nearby tracks. We 
therefore define 

£(c£) = • £mc(c£), (58) 

where £ norm is the normalization needed to scale the simulation to agree with the data. 

The topology of a B — > u£D** decay chain is shown in Fig. |S[. The ct dependence of 
£ is the result of the impact parameter significance cut (d /o~ < 3.0) in the SST selection 
(Sec. [IV B| ). By removing this cut, we remove the ct dependence from £. Figure [Tt] shows 
£,Mc{ c ~t) with the d /a cut removed (top), and with the cut applied as normal (bottom). 
Without the cut the distribution is flat, as expected, and corresponds to a 33% probability 
to tag on a nf^ given that one is present. Applying the do/&o cut rejects most of the 7r^ tags, 
especially once ct is beyond a few hundred microns. The £ shape is modeled by a Gaussian, 
centered at zero, with a constant term. 

To determine the normalization, £ norm , we remove the d /a cut from the data (analogous 
to Fig. [K| top), thereby eliminating the ct dependence as well as enriching the sample in 
T^t* tags. We divide the tagged events into right-sign and wrong-sign tags, and make the 
distribution of the impact parameter significances with respect to the B vertex {dsjoB of 
Fig. H). An example of such a distribution is shown in Fig. [TTJ. The excess of right-sign events 
near d^jo^ = is due to the 7r^ tags. Their number, iV(7r*" 9 ), is determined by fitting the 
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distribution to a Gaussian (centered at zero and with a unit RMS) for the 7r^'s, plus a nf^ 
background shape obtained from the wrong-sign distribution. The wrong-sign distribution 
renormalized to the fit result is overlaid onto the right-sign distribution in Fig. O. It is 
seen to agree very well with the right-sign distribution at large de/crg, and displays a clear 
excess near zero. 

We also count, again without the c?o/co cut, the total number of tags N(tags) and 
compute the fraction of B decays where vr^'s are tags, 

N(ir tag ) 

Rk = N^j- (59) 

for signature k. The measured ratios are given in Table 0. We extend our notation so that 
%mc is £mc( c ^) when the impact significance cut, d /<j < 3.0, is removed. %mc is then 
independent of ct. The _R**'s are simply 

POO 

Rk = UrmLlC / {If® + T^\t)} dt, (60) 
JO 

for decay signature k. Thus, £ n0 rm is simply related to R^*, the other sample composition pa- 
rameters, and %mc- Rather than attempting to compute an average £ n0 rm, we will constrain 
the five R* k * predictions to the measured ratios in the Arrid fit (Sec. |V F 1| ) . 



F. Fitting the asymmetries 

1. The x 2 function 

The observed tagging asymmetries can be predicted in terms of the sample composition 
parameters and the true asymmetries. The true asymmetry for the B + is constant in ct, 
while for B° it follows a cosine dependence, and accounting for the ct resolution one has 



A + (ct) =V + (61) 
o _ G{ct; ct', a ct ) g {e-t'/^Vp cos{Am d ct')} 



where ® denotes convolution of the physical time dependence (cosine and/or exponential 
functions) with the ct resolution function G(ct; ct', a ct ) over ct'. The latter function is a nor- 
malized Gaussian of mean ct' and RMS a ct . However, the measured ct, and associated reso- 
lution, depends upon the sample composition. Therefore, the proper times used for the pre- 
dicted asymmetries are the ct° k ,+ [Eq. (|49|) 1 obtained using the sample composition-averaged 
/C-factor. For the resolution a ct we use the composition- weighted resolution a k (ct Q k ,+ ) from 
Eq. ©. 

We form a \ 2 function to simultaneously fit Am^, T> + , and T>o over all ct-bins of all 
decay signatures by comparing the predictions Ak{ct kd ) calculated via Eq. (|5T|) against the 
measured asymmetries A k meas \ct kd ) , where ct kd is used for j\t™ eas ^ since we were restricted 



to the direct chain ctkd when binning the data (Sec. |VD|) . The Ak(ctkd) asymmetry depends 



not only upon the parameters Am^, Vq, and T> + , which are of direct interest, but also on 
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t , t + /t , /**, Py, e(7r*), and ^ nor m through the JF° i+ 's. The last two parameters are also 
expressed as functions of R* and R%*, as well as the other composition parameters. The 
comparison of A k meas ^ and Ak corresponds to the first summation in our x 2 function: 

A k meas \ct k d) - A k (ct kd )^ 




F ^ -F 3 {r,P v ,R f ,^ s )...) 
-I 



(63) 



where k is an index that runs over the five decay signatures, and ctkd symbolizes the sum- 
mation over the proper time bins. 

The second summation is over the set of fit input parameters: p\ meas ) j g 
measured input value for parameter j, a J is its error, and the "predicted" value is 
Fj(f**, Py, Rf, e(7r*), . . .). This prediction is a function of the sample composition parame- 
ters, and in most cases it is a trivial substitution, such as Fj(f**, Py, Rf, e(7i*), . . .) = r for 
the B° lifetime. However, e(7r#), (norm, and Py are not directly measured but are constrained 
in the fit by their appearance in the predictions for other measured quantities, namely R* 
and R* k * . Allowing e(7T«), (.norm, and Py to float in the fit constrained by the measured 
sample composition parameters was one of the motivations for extending the x 2 function 
with the second summation. 



The reconstruction efficiency for the D* pion e(-7r*) can be obtained (Sec. |VB 3|) from 
R* ( meas \ measured in the data to be 0.249 ± 0.008. Since the prediction R* is a function of 
the sample composition parameters, e(7r*) depends on them also [Eqs. (f43"l)-(f45l)1 and must be 



recomputed whenever the sample composition parameters change. This recomputation nat- 
urally occurs in the \ 2 minimization by allowing the composition parameters that determine 
R* to float, coupled with the constraint of the term 

R* {meas) -W(e(^),r,...) ^ 2 (64) 

in the x 2 . 

A similar approach holds for (norm and Py using the R**'s. In this case there are five 
.RJf's, one for each decay signature, and a Fj term for each. The prediction for R* k * is 
proportional to (norm by Eq. flBT^). Because (norm is common across decay signatures, it 
is essentially determined by the average of all five R** {meash s. Py, the relative D** decay 
rate to D* vs. D [Eq. (|35|)1, is treated as completely unknown. However, it is also related 
to the -Rfc*'s. If Py = 0, there would be no D** — > D* decays, and consequently no 7r„'s 
in the £ + D*~ signatures, resulting in the corresponding R* k = 0. The values of the -R£*'s 
relative to each other determine Py. While the errors on R* k * (meas ) are large (Table 0), and 
therefore Py is not tightly constrained, this method is preferable to just using a theoretical 
estimate for Py. Its incorporation into the x 2 function automatically enables it — like the 
other parameters — to vary within the allowed experimental constraints and propagates the 
associated uncertainty to the fit parameters. 
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2. Result of the fit 



The x 2 function is minimized over the six ct bins for all five decay signatures simulta- 
neously, letting the unknown parameters float freely, and the known inputs to vary within 
their errors. The fit results in the following values 

Am d = 0.47li°'°?5 P^ 1 
V + = 0.267j£gS 
X?o = 0.181^6 

efa) = 0.845lffi 
p = 747+ - 470 

S norm u . i *± i _ 



-0.292 

Py = 0.33ltollt, 

with the nominal fit errors quoted. The x 2 °f the fit is 26.5 for 30 degrees of freedom. 

The 7T* efficiency is quite high, thereby limiting one source of the B° <-> 5 + cross-talk. 
The D** contribution to cross-talk is quantified by Py, which is on the low side of what is 
sometimes assumed |26| . Our value could be biased by our sample selection, but in any case 
the errors are large. The final sample composition yields: ~ 82% of the £ + D°X signature 
comes from B + decays, while ~ 80% of the £ + D~X and ~ 95% of the £ + D*~X originate 
from B°. 

Figure O shows the result of the fit overlaid on top of the measured asymmetries, where 



all three £ + D*~ signatures are combined. Figure [13] gives the three £ + D*~ signatures sep- 
arately. The cross-talk is relatively modest, and the signatures dominated by B° gener- 
ally show a fairly clear oscillatory behavior in the raw observed asymmetries. For the 
i? + -dominated signature, £ + D°, the raw asymmetry is compatible with being a constant 
(Fig. [12p, but the residual effect of the B° cross talk is visible in the fitted curve in the form 
of a slight oscillation. The effect of the B + contamination in the B° signatures amounts to 
a constant shift in the asymmetry and is therefore not readily apparent. 



The fit parameters constrained to a priori measured values are shown in Table [Vj along 
with the value and error output by the fit. Except for /**, these parameters are largely 
decoupled from the other fit parameters, and are virtually unchanged from their input 
values. The data are more sensitive to the value of /** because it governs the amount of 
cross-talk between B° and B + . 

The correlation coefficients of the fit parameters with Am^, T> + , and V are shown in 
Table |V11| . We see that the lifetimes are largely decoupled from other parameters. On the 
other hand, /**, (, n0 rm, and Py are strongly coupled to Arrid, T> + , and T> , underscoring the 
importance of the 7r^; corrections to the analysis. The correlation between Am^ and T> + is 
stronger than between Am^ and D . The reason for this stronger corrrelation is that the 
effect of T> + enters via the B + contamination of the B° signatures and is manifested by a 
downward translation of the B° oscillation in Fig. 12]. As the oscillation is translated down, 
the intercept with zero asymmetry moves to shorter times, thereby decreasing Am^. On the 
other hand, if Dq is varied, the oscillation amplitude varies symmetrically about the vertical 
axis and Am d is weakly affected. 
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G. Statistical and systematic uncertainties 



We fit for Amj, T> + , and T> using a y 2 function which also incorporates the sample 
composition parameters. The errors it returns are a combination of statistical and systematic 
effects, yet the errors only partially account for the systematic uncertainties. The sources of 
the uncertainty can be divided into statistical and three systematic categories: 

• Statistical: the error that is directly correlated with the £D^ sample size. 

• Correlated Systematics: parameters of the fit (/**, Rf, t + /t and r ), coupled to 
Arrid, T> + , and V through the sample composition [Eq. (f)7|)]. These parameters are 
not correlated among themselves; only their effects on Am^ T> + , and V are. 

• Uncorrelated Systematics: systematic uncertainties caused by imperfect simulation 
models of the physics processes or the detector. 

• Systematics from Physics Backgrounds: uncertainties due to other physics pro- 
cesses that contribute to B — > u£D^*' data samples that have been hitherto neglected. 

We determine the uncertainties from each of these four categories in turn to estimate the 
statistical and systematic uncertainty for our final result. 

We separate the statistical and correlated systematic errors of the original fit by repeating 
the fit with the sample composition parameters fixed to the results of the original fit, and 
only six variables (Am d , V + , V , e(7r s ), („ orm , and Py) floating. The errors from the six- 
variable fit are just statistical (cr stat ), while the errors from the full fit are the combined 
statistical and correlated systematic errors (cr s tat+c.s.) ■ hi a Gaussian approximation, we 
estimate the correlated systematic error, uc.s., by 



<?C.S. = \J Vstat+C.S. - Vstat I (65) 

and find, for example, 

A™ _ n 471 +0.078 +0.030 -1 
am d _ u -^' i -0.068 -0.031 P b J 

where the first error is purely statistical and the second is the correlated systematic. This 
correlated error is listed in Table |V111| under "Sample Composition," and is by far the 



dominant systematic uncertainty. 

The "uncorrelated" uncertainties include the contributions from the uncertainty in the 



Monte Carlo modeling, which are also listed in Table |VIII| . An uncertainty in the 6-quark 
production spectrum (App. |A 1|) carries over into the determination of the /C-factor dis- 
tributions. The systematic uncertainty was estimated by weighting the generated pr(B) 
distribution by a power law factor whose range was obtained from a b cross section analy- 
sis using an inclusive electron sample. The shifts in fit parameters using this weighted pt 
spectrum are taken as the associated uncertainty. 

The isolation requirement of the inclusive electron trigger (i.e., no matching cluster 
in the hadronic calorimeter), could, if poorly simulated, bias the decay kinematics of the 
selected B's, and result in an erroneous /C-factor. Since this requirement is not present in 
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the inclusive muon trigger, we use the difference obtained in the fit when using the electron 
sample composition parameters versus those of the muon for this uncertainty. 

Various calculations (e.g., of efficiencies, /C-factors. . . ) are sensitive to the differences in 
the B decay dynamics to a D, D* , or D** . The systematic uncertainty due to the decay 
model is obtained by repeating the analysis where the decays are governed by phase space 
instead. 

The L xy resolution used in this analysis is from the CDF detector simulation. The 
systematic uncertainty is obtained by varying the intrinsic resolution by ±20%, and the 
resultant shifts are taken as the uncertainty. 

The last uncorrected systematic error is due to the shape of £a/c( c ^)? the time depen- 
dence of the probability to tag on a n„ from D** decay. An alternative shape for £mc*(c£) 
is obtained by using the £D*~ signature instead of £D~ , and using a variant of the CDF 
detector simulation. The £mc( c ^) is again well described by a Gaussian plus a constant, but 
the new RMS of the Gaussian is 400 /im, which is twice the nominal value. The shifts in the 
fit resulting from this wider ^Mc{ct) are taken as the systematic uncertainty. 

Our results may also be affected by physics backgrounds not included in the sample 
composition which result in uiD^X with the correct correlation of £ and D^*': 

• B -> DPDWX, followed by -> v£X; 

• B s -> v£D**, followed by D** -> D^K; 

• gluon splitting g — > cc, followed by c — > £X and c — > D^Y. 

The fractional contributions of the first two processes to our sample are estimated by 
Monte Carlo simulation (App. |A f| ). The fractions, listed in Table fX|, are small. 

Because of the uncertainty in accurately predicting the rate and other characteristics of 
gluon splitting, we use data rather than simulation to set an upper bound on this contribu- 
tion. For this background the apparent B vertex is reconstructed from two different charm 
decays, so the reconstructed D will have a broad etc distribution, including cases where the 
D apparently decayed before the U B." There is some (statistically marginal) evidence of 
right sign (^K^) signal in the cto < —1 mm region in the data. We use the size of the far 



negative cto tail to constrain the potential size of the gluon contribution |27_]. Because of 
the large statistical uncertainty, we take as the upper bound on the gluon contribution the 
central value of our fitted fraction plus twice the statistical error on the fraction (Table |TX| ). 
Doubling the statistical error is ad hoc, but we wished to be conservative in accounting for 
this poorly constrained process. 

The effect of each physics process on the asymmetry is accounted for by adding two 
new terms to the predicted asymmetry of Eq. ([57]), one for tagging on fragmentation tracks 
A^ ra9 , and another for tagging on decay products A decay . We can determine each of these 
asymmetries, or their upper bounds, and combined with the composition fractions repeat 



the fits to the observed asymmetry |27j. The shift in fit output under each of these processes 



is taken as their contribution to the systematic uncertainty. 



Examination of Table |VIII| shows that by far the largest contribution to the systematic 
uncertainty comes from the input sample composition parameters. The combined systematic 
uncertainty is obtained by adding the individual contributions in quadrature. The combined 
systematics are still smaller than the statistical uncertainties, especially in the case of Am^. 
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As a mixing measurement, the application of Same Side Tagging on the £D^*> sample is still 
limited by statistics. 



H. Summary of the ID^*> analysis 

We have applied our SST method to a partially reconstructed B — > £D^ sample and 
accounted for the effects of B° «-> B + cross talk in the sample composition. The flavor 
oscillation is readily apparent, and the oscillation frequency and dilutions are found to be 



Am d = 0.471 

P4 



hO.078 
-0.068 



± 0.034 ps" 



0.267 ±0.032 tomi 
V = 0.181 ±0.028 toml 



(66) 
(67) 
(68) 



where the first error is statistical and the second is the combined systematic 
value compares well with a recent world average of 0.484 ± 0.026 ps _1 |29] 



Our Amd 
We will discuss 



the dilutions further in Sec. VII 



VI. TESTING SST IN B -> J/ipK SAMPLES 

Having demonstrated that our SST algorithm is capable of revealing the B°-B° fla- 
vor oscillation in a large lepton±charm sample, we extend its use to the exclusive modes, 
B — > J/ipK + and J/ip K*° (892), where one may further test this method. The SST dilu- 
tions should, ignoring some experimental biases, be independent of decay mode, and the 
J/ipK samples should yield results comparable to the £D^*> analysis. These J/ipK samples 
are too small to provide precise tagging measurements, but they provide an experimental 
opportunity to study flavor tagging in this type of exclusive mode. This study is especially 
interesting because it serves as a model for tagging B° — > J/ipK^, which we consider in 
Ref. [Q. 



A. Reconstruction and tagging of B — > J/ipK 



Our J/i/jK samples have appeared in a number of previous publications, in whole or part, 
on measurements of B masses |[30|| , lifetimes ]T5| , |3T| , branching ratios and production 
cross sections ||33||. The reconstruction criteria are somewhat different here; we wished to 



maximize the effective statistics and were less concerned about accurately modeling efficien- 
cies or triggers. 

B reconstruction begins by forming charged particle combinations with J/ip candidates 
(Sec. |I1 C| ). Since we require a well measured B vertex, at least two particles of the decay 
must be reconstructed in the SVX with loose quality requirements (principally that the track 
used hits on at least 3 out of 4 SVX layers). For the J/ipK + a single particle, assumed to 
be a kaon, with px > 1.75 GeV/c is combined with the J/ip. The J/ipK*° reconstruction 
uses pairs of oppositely charged particles, each with p T > 0.5 GeV/c. The pair is accepted 
as a K*° candidate if a vertex-constrained x 2 fit — considering both permutations of K ± 
and ir T mass assignments — yields a mass within 80 MeV/c 2 of the K*° mass (896 MeV/c 2 ), 
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and has pt{K*°) > 3.0 GeV/c. The fit includes dE/dx energy loss corrections appropriate 
for the mass assignments. If both permutations satisfy these requirements, the assignment 
closest to the K*° mass is selected. The high pt(K*°) cut is necessary to reduce the larger 
combinatoric background for J/ipK*°. 

The particle(s) making the K + (K*°) are combined with the pair in a multiparticle 

X 2 fit for the B with the mass constrained to the world average J / ip mass, all daughter 

particles originating from a common vertex, and the entire system constrained to point to 
the pp interaction vertex. A run-averaged interaction vertex is used as was done for the 
£D^ sample (Sec. |V A). The RMS spread of the transverse beam profile is taken to be 40 



/im. The resulting B candidate must have Pt(B) > 4.5 GeV/c. 

Rather than cutting on the \ 2 from the multiparticle fit, we cut on only the portion 
coming from the transverse (r-0) track parameters {i.e., curvature, azimuthal angle, and 
impact parameter The B pointing resolution to the primary vertex in the r-z plane is 

very coarse in CDF, providing little separation between signal and non-pointing backgrounds. 
Including the z-cot 9 contributions to the \ 2 tends to smear the separation that is available 
from the precise transverse measurements. We require the transverse tracking terms of the 
X 2 sum to be less than 20, and similarly that the transverse components of the vertex \ 2 
sum to be less than 4. Although formally ad hoc, we found these cuts to be a little better at 
discriminating signal from background than the full y 2 . However, the final tagging analysis 
is insensitive to the type and value of the x 2 cu ^s used in the B reconstruction. 

Finally, if there are multiple B candidates in the same event, the one with the smaller 
transverse track parameter \ 2 is taken. 

These B samples are used in a likelihood fit (Sec. |V1B| ) employing a normalized mass 
variable M^, and so we discuss the selection results in terms of this variable. is defined 
as (MpiT — M )/<t fit , where M FIT is the mass of the B candidate from the fit described 
above, M is the central value of the B mass peak (5.277 GeV/c 2 ),0 and a fit is the mass 
error from the fit. Over the range of | |< 20 we have a total, signal plus background, of 
12564 events in the J/ipK + sample and 2339 for the J/tpK*°. 

Figure [14] shows the normalized mass distributions for candidates with reconstructed 
ct > [Eq. ([46"D1. Also shown is the result of the likelihood fit performed in Sec. |V1D| , where 
the mass is modeled by a Gaussian signal plus linear background. The fit yields 846 ± 35 
S+'s and 365 ± 22 B 0, s (for all ct). 

Events with ct < are dominated by background (see Fig. [16| in Sec. |VI D| ), and the 



mass distributions show no clear signal. However, these events help constrain the background 
behavior and are kept as part of the analysis. 

These two samples are then tagged with the SST criteria of Sec. [IV H . We find about 
63% of the J/ipK + and J/ipK*° signal events are tagged. 



8 The mean is systematically low by 2 MeV/c 2 compared with the world average mass because we 
do not make all the detailed corrections of Ref . fpOl . 
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B. Likelihood function for the J/ipK samples 



The tagging correlations of the B — ► J/ipK samples have the same physical time de- 
pendence as the £D^ samples [Eqs. and (§)] but without the complications of sample 
composition and average /37-corrections. Maximal use of the smaller J/ipK samples moti- 
vated a more sophisticated approach than used to fit the £D^ data. 

An unbinned likelihood function is devised to simultaneously fit over various measured 
event properties and obtain the SST dilutions for the J/ipK samples. The likelihood function 
incorporates the candidate's proper decay time and invariant mass to facilitate separation of 
signal and background. It is also generalized to consider tagging biases. These are relatively 
unimportant in mixing measurements which only use the relative flavor-charge asymmetry 
but are critical for CP violation measurements where the effect appears as an absolute 
charge asymmetry of the tag. Although our focus is on the charge-flavor correlations of 
SST, this generalized approach serves as a prototype for CP violation studies 

The likelihood function to be maximized is given by 



N 



£ = U^ bCb + - fB)(hC L + (1 - fi)£p)} (69) 

8=1 

where the product is over all N events in the mass window |Mjy| < 20. The subscripts 
B, L, and P respectively indicate terms associated with the 5-meson signal, long-lived 
backgrounds, and prompt backgrounds. The fraction of events that are B signal is fg. 
The backgrounds are subdivided into two classes: "long-lived," which are those consistent 
with a non-zero lifetime, and "prompt," which are those consistent with zero lifetime. The 
fraction of long-lived backgrounds, which are predominantely real 5's that have been mis- 
reconstructed, is given by 

The (<f) = B, P, and L) are functions describing the relative probability for obtaining 
the following measured values in an event: the normalized mass (Mat), the proper decay 
time and its uncertainty (t and at), the reconstructed decay flavor r (r is +1 for B + and 
B°, and —1 for B~ and B°), and the tag track sign s (s is +1 for a positive track, —1 for a 
negative track, and if there was no tag). 

The density function for the J/ipK + signal describing the mass and t dependence, the 
relative numbers of B + and B~, the flavor-charge tagging correlation, and the tagging 
efficiency is 

C B = G(M N ;0,X) G(t;t',Ya t ) ® E(t';r B ) X 

l + rR B \ (1-™ B {s)V b \ Eb[s) {7Q) 

G(x; [A, a) is a normalized Gaussian distribution in x G (—00, +00) with mean ll and RMS a, 
and E(x; t) is a normalized exponential distribution in x £ [0, +00) with mean r. The first 
factor in C B is the shape of the mass distribution (Mat), where A is a scale factor for the 
mass error. The second factor is the Gaussian resolution of the reconstructed t, including a 
resolution scale factor Y . The <g> denotes convolution over t', in this case with an exponential 
distribution E of lifetime r B = r+. The resolution scales, A and Y, are extra degrees of 
freedom to monitor our description of the errors. 
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The density function next contains two asymmetry factors. The first is the probability 
of reconstructing the observed meson flavor r, and depends upon 



Rb 



N(B+) - N(B~) 
N(B+) + N(B-)' 



(71) 



This first factor decouples other flavor-related asymmetries from a "reconstruction asym- 
metry" Rb by accounting for the different numbers (N) of 5 +, s and B~ J s that may be 
reconstructed due to a detector bias, or simply a statistical fluctuation in the relative yield. 

The second asymmetry factor represents the probability of obtaining a tag of sign s given 
the reconstructed flavor r. The strength of the flavor-charge (r-s) correlation is the usual 
dilution, T>b = T> + . The effective tag for a track of sign s is k_b(s). Since the B-tt correlation 
is between B (r = ±1) and tt t (s = =Fl), the tkb{s) term appears with a negative sign. 
Finally, the efficiency to obtain such a tag is £b{s). 

This formulation with T>b, kb{s), and £b{s) is able to account for the general situation 
where the tagging method suffers from intrinsic tagging asymmetries, as might be caused by 
detector biases. Tagging asymmetries may result in different dilutions and efficiencies for the 
two B flavors. We define V>b to be the flavor- averaged dilution, and we are able incorporate 
all tagging asymmetry effects in Kb{s) and £b(s). In the absence of tagging biases, «b(s) is 
simply the charge of the tagging track (kb(s) = s), and £b(+1) = 1) = e B , where e# 
is the flavor-averaged tagging efficiency. 

The "charge asymmetry corrected" tag kb{s) and efficiency £#(s) are determined using 
two new parameters: as, which is the charge asymmetry in selecting a tag track, i.e., a 
bias in selecting positive vs. negative tags; and 5b, which is the 6-flavor asymmetry in 
finding a tag track, i.e., having different efficiencies to tag b vs. b mesons. For convenience 
we sometimes normalize the latter by the dilution, 75 = 8b/ Db- The derivation of kb{s) 
and Sb{s) may be found in App. [B|, along with a complete characterization of a generic 
tagging method. The actual determination of the tagging biases in our detector is discussed 
in Sec. gTUg 

The density function for the J/ipK*° signal takes the same basic form as for J/ipK + , 
but it now incorporates the time- dependent flavor oscillation A°(t') = cos(Amdt'). Since 
no particle identification is used, some J/ipK* events enter the sample with the correct 
K-tt mass assignments "swapped" (Sec. [VI C|) , for which the apparent flavor is inverted. 
The density function is divided into unswapped and swapped parts, with the reconstructed 
flavor r for the swapped events appearing with a reversed sign. The complete expression is 



where f$ is the fraction of swapped events, and /x$ and X$ are the mean and RMS of the 
normalized mass distribution for the swapped events. The rest of the parameters parallel 
those of the J/ipK + , but with V B = T> . 



Cb 



(1 - f s )G(M N ;0,X) G(t;t',Ya t ) ® [E(t';r B ) x 
(l+rR B \ ( l+rKB(s)V B Mt')\ ' 




+ f s G(M N ; t x s ,Xs)G(t;t',Ya t ) ® [E{t';r B ) x 
(l-rR B \ f 1 - rK B (s)V B Ao(t')\ ' 




(72) 
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The density function for the long-lived background for both decay modes is similar to 
the signal except for a linear mass distribution, the presence of three exponential lifetime 
distributions, and the lack of mixing: 



f l + ( L M N \ 

{ 2W ) G fr*' y<T «) 



{f N E(-t'; r L2 ) + (1 - f N )[f T2 E(t'; r L2 ) + 
(l-/ T2 W;r L i)]}x 
(l + rR L \ fl + rK L {s)V L \ 




(73) 



The linear mass distribution is parameterized by a slope and the width of the mass 
window | | < W = 20. The long-lived background consists of positive- and negative-t 
components, with a fraction f N in the negative exponential (with lifetime r L2 ). The positive- 
t background is described by two exponentials, one with a large lifetime t l1 , plus a short 
one of Ti2- The latter lifetime is fixed to be the same as for the negative-t tail. The fraction 
of positive-t events (1 — /jv) which compose the short lifetime exponential is f T2 - 

The background may also possess a reconstruction asymmetry Rl, or a charge correlation 
between the tag and what is assumed to be the K + or K*°, i.e., a dilution V L . The 
background asymmetry description parallels that of the signal with R^, D^, K<p, and 
defined independently for each event class (0 = B, L, and P). 

The prompt background density function is 



with the same variable definitions as before except that they apply to the prompt back- 
ground. There is only a dependence on the proper time through the t resolution, and thus 
no convolution is needed. 

When summed together and multiplied over all the selected events in a particular dataset, 
the density functions form a properly normalized likelihood function. 



A number of the likelihood parameters are more accurately obtained from sources other 
than our J/ipK data. In this section we will discuss which parameters are fixed in the fit, 
and their sources, values, and uncertainties. 



The likelihood function relies on the temporal properties of the B decay, and these are 
best obtained from world averages. Since we wish to measure the tagging dilution, and not 




(74) 



C. Input likelihood parameters 



1. B meson parameters 
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the oscillation frequency, we include Am^ in this list. We use the following averages from 



the Particle Data Group |22 



Am d = 0.474 ± 0.031 ps" 1 
r+ = 1.62 ±0.06 ps 
t = 1.56 ± 0.06 ps. 



2. Incorrect K — ir assignment 

The J/ipK*° events include real B° —>■ J/i/jK* decays, but with the incorrect K-n 
mass assignment. A Monte Carlo sample of B° decays (App. [A 1| ) was generated and then 
processed as data. The reconstruction tries both assignments, and if both versions pass 
the selection criteria, the one with the Kit mass nearest the K* 0, s is chosen. The events 
with the K and 7r swapped have an Mjy distribution which is roughly Gaussian with mean 
Us = —0.5 and RMS Xs = 4.9. The area of the swapped Gaussian is 9.8% of that for the 
unswapped distribution. The kinematic dependence of the swapped events on Pt(B) has 
also been studied ||34j| . The fraction of swapped events is constant within a few percent over 
our range of pr{B), but the mean and RMS of Mjv show some systematic variation. 

The swapped component is difficult to fit in the J/i/jK*° data because it is difficult to 
distinguish these events from the combined shapes of the narrow central Gaussian and linear 
background. The likelihood fit therefore fixes the input parameters /15, Xs, and fs to the 
values from the simulation. We allow for a 100% variation in the fraction, fs = 0.1 ± 0.1, 
and assign uncertainties to the other swapping parameters which covers the range observed 
when spanning the Pt{B) interval of the data, i.e., fis = —0.5 ± 0.5 and Xs = 5.0 ± 2.0. 



3. Tagging biases 

A tagging method may display two sorts of inherent asymmetries (App. 0): selecting 
one charge more often than the other as a tag (a), or having a greater efficiency to tag on 
one 6-flavor over the other (<5). 

The charge asymmetry of the tags for a flavor symmetric sample is 

N+-N-. 

a = — — , (75) 

iV + + iV_ v ; 

where N + (iV_) is the number of positive (negative) tags. Since we reconstruct the decay 
flavor, we can correct for any flavor asymmetry in our samples and determine a from the 
data. This is done for the backgrounds by letting and ap float in the likelihood fit. 

However, a appears in the likelihood function partly as a factor T>/(1 ± a) [Eqs. ( |70|) 
and (PTD1- We can essentially eliminate the influence of the T>b-o>b coupling in the fitted 
T>b by fixing as to an independently measured value; we obtain a better constraint on a# 
as well. We do this by using a large inclusive sample of non-prompt J/^'s, i.e., a flavor- 
symmetric b sample. This is the sample of J/^'s described in Sec. |II C| with the following 
additional requirements: both muons are in the SVX, and the J /if) projected flight distance 
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Lxy(J/if>) [Eq. (|T^) 1 exceeds 200 fim. This last cut results in a sample which is more than 
90% pure b hadron decays. We also require pr(J/if)) > 4 GeV/c so that the pts are similar 
to that of the J/ipK samples. 

We have looked for tagging asymmetry dependencies in a variety of variables, including 
the pt and L xy of the J/ip, and found only two variables of interest. First, there is an 
a-dependence on the pt of the tagging track, with more positive tracks reconstructed than 
negative ones at low pt- This is due to the charge asymmetry inherent in the design of the 
CTC (wire planes are oriented along the direction of positive tracks). Proton spallation from 
the beampipe might contribute an additional p^-dependent asymmetry, but this effect has 
been largely eliminated by the impact parameter significance cut on the SST candidates. 
The second variable is the number of good pp interaction vertices ny found by the VTX. The 
number of vertices is an indicator of the total hit occupancy in the CTC, which influences 
the tracking efficiency. 

Characteristics of the J/ ipK samples, or the criteria of the tagging method, could modify 
the biases in the tracking asymmetry. We therefore compare the charge asymmetry in four 
types of J/ip (L xy > 200 /xm) subsamples (the J/ip is used for the b direction): (i) the SST 
tags, (it) SST candidates (i.e., no pTp 1 selection), (tit) SST candidates passing a 6- vert ex 
veto, and (iv) all tracks satisfying the SST cuts except for being in a "side cone" away 
from the J/^-axis. The first case is the most direct extension of our analysis, but it is 
contaminated by tagging on B daughters. The second case has greater statistics as there 
are multiple track entries per J/ip, but more importantly there is no bias associated with 
the fir 1 cut. The third sample suppresses tagging on B daughters — such tags are impossible 
with the exclusive J/ipK signal — by requiring that the SST candidate impact parameter 
significance relative to the J/ip vertex, d^/a^, is greater than 2. In the final case we select 
tracks with the basic SST cuts, but which are divorced from the b by using all tracks in a 
1.0 < AR < 2.1 "side cone" relative to the J/ip. 

Since we wish to study a pr dependence, and only sample (t) selects a unique track per 
J/ip, we relax the 400 MeV/c SST cut on samples (ii)-(iv). The upper plot in Fig. [L5| shows 
the charge asymmetry as a function of the track's inverse pt for these four samples. The 
asymmetry is fairly small at 1/pr = 2.5 (GeV/c) -1 (the nominal SST p? cut of 400 MeV/c) 
but then rises significantly. The asymmetry as a function of number of vertices ny is shown 
in the lower plot in Fig. |Tj| Neither variable shows any significant differences across the 
four samples. 

With all four samples so similar, we consider the asymmetry to be independent of the 
tagging, and choose sample (Hi) to determine a parameterization for the tag asymmetry a^. 
The pt dependence of the asymmetry is well described by p^, and a linear function for the 
number of vertices. In terms of both variables we write 

®b(Pt, n v ) = {a x (ny - n ) + b 1 }(p7 r A - p~*) 

+ {a 2 (n v -n ) + b 2 }, (76) 

where we have included offsets p^ and n , fixed to 20.0 (GeV/c)~ 4 and 3.0 respectively, to 
remove the correlation between the a's and 6's when fitting for them. 

We determine these parameters by making subsamples for each integer value of ny, and 
fitting them for the coefficient of the p^ 4 — p^ term and a constant offset. The series of 
these coefficients and offsets are then fit for the linear ny dependence. We obtain the values 
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ai = (3.9 ± 1.8) x 1(T 4 (GeV/c) 4 
h = (1.3 ± 0.4) x 10" 3 (GeV/c) 4 
a 2 = (1.4 ±0.4) x 10~ 2 
b 2 = (2.6 ±0.8) x 10" 2 , 

which give the curves in Fig. |15[ If we consider the tags in the J/ip data [case (i)}, we 
find the average tag asymmetry is (1.6 ± 0.7)%. In the limit as 1/pr - » the asymmetry 
parameterization gives (0.14 ±0.86)%. We use this parameterization to describe the tagging 
asymmetry for the B signal as in the likelihood fit (Sec. |VI D| ). 

As well as an intrinsic bias towards positive or negative tags a, the SST could also have 
a bias 5 to tag b and b mesons with different efficiencies (see App. [FJ), i.e., 

S = 7mT7m- (77) 

where e is the efficiency to tag on a given flavor. It is more convenient to express its ratio 
relative to the dilution in the likelihood, so we often use 7 = 5/T>. We can constrain 5 from 
the data, and do so for the backgrounds by letting 5p^ float in the likelihood. 

In the likelihood function 7^ appears partly as a factor (1 ± j B )V B . As with a B , an 
independent determination of 75 is preferable in order to decouple it from the dilution. 
The inclusive J/ip's cannot be used since we have no knowledge of the b flavor. The higher 
statistics £D^ data indicate that jb is less than 15%. However, as discussed in Appendix 0, 
we can improve upon this constraint by about a factor of three by considering the behavior 
of •ys/ciB- We find ^bI&b spans the range from 0.0 to about 2.5, and we use •Jb/c^b = l-Ol^o 
in the likelihood fit. 

This completes the list of input parameters fixed in the likelihood fit, and we now proceed 
to fitting the data. 



D. Fitting the J/ipK + and J/il)K* Q samples 

We use the likelihood function to fit the J/%l)K samples with the parameters discussed 
in the last section fixed, and allow the others to float freely, i.e., 

• fs, the fraction of events which are signal, 

• fi, the fraction of background which is long-lived, 

• /iv, the fraction of the long-lived background in the negative lifetime tail, 

• f T 2, the fraction of positive long-lived background with lifetime tl2, 

• Cp and (l, the slopes of the prompt and long-lived backgrounds in normalized mass, 

• X, the error-scale factor for the normalized mass, 

• Y, the error-scale factor for the decay time, 

• r L1 , the large lifetime of the positive long-lived background, 
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• tl2, the small lifetime used for positive and negative long-lived backgrounds, 

• €b, €p, and el, the efficiencies for tagging signal, and prompt and long-lived back- 
grounds, 

• Rb-i Rp-i and Rl, the reconstruction asymmetries of the signal, and prompt and long- 
lived backgrounds, 

• atp and «l, the tagging asymmetries of prompt and long-lived backgrounds, 

• 5p and Sl, the tagging efficiency asymmetries of prompt and long-lived backgrounds, 
and 

• Vb, T>p, and Vl, the dilutions for the signal, and the prompt and long-lived back- 
grounds. 

Each decay mode is separately fit by minimizing the negative log-likelihood function. 

The fit results are compared to the distribution in Fig. |T3|. The fitted proper decay 



lengths are shown in Fig. [16] for the J/ipK*°, where we have defined — for display purposes 
only — the signal region as |M/v| < 3, and 3 < \M^\ < 20 as sidebands. The data are well 
described by the fits. 

To display the flavor-charge asymmetries, we compute the mass sideband subtracted 
asymmetry (analogous to Eq. (|55|)) for the data in ct bins. The results are shown in Fig. |H] 
with the likelihood fits superimposed (solid line). The J/ipK + plot shows a clear correlation, 
consistent with being constant, and the J/ipK*° data is in good agreement with the mixing 
hypothesis. Also shown in the figure insets are —2 hi(£/ £ MAX ) as a function of the dilution, 
where £ is the value of the likelihood for a given T>p after maximization with respect to all 
other free parameters, and Cm ax is its value at the global maximum. We see a well behaved, 
approximately parabolic, shape. The 1, 2, and 3a errors of the likelihood are indicated in 
the inset by the three horizontal dotted lines. 



As a simple check, the binned and mass sideband subtracted asymmetries of Fig. [T7| were 
fit to a constant for J ji\)K + , and T>q cos(Amrft) for J/ipK* , using a simple x 2 fit without any 
additional corrections (e.g., tagging asymmetries, t resolution, etc.). The results (dashed 
line in Fig. |T7|) agree very well with the likelihood fits, indicating that the fits are driven by 
the basic asymmetries in the data and are not significantly influenced by the refinements of 
the likelihood fit. Of course, the fits are strongly dominated by the statistics. 

The principal results of the likelihood fits are the dilutions T>b given in Table |X|, along 
with the other fitted parameters. A few remarks may be made on this table in passing. The 
t error scales Y are virtually unity, indicating that the lifetime modeling and error estimates 
describe the data well. The X scales (the RMS of the signal M N Gaussian) are not 1.0, 
but are instead close to the known scale of 1.3 ||30|| . The tagging asymmetries a for the 
background terms are generally consistent with the (1.6 ± 0.7)% found as the average value 
from inclusive J/^'s (Sec. |V1 C 3| ), with otp for the B 0, s about 2.5a high. The reconstruction 
asymmetries R are also not statistically significant beyond l-2er. The background dilutions 
are, not surprisingly, consistent with zero when selecting pairs of tracks with no net charge 
(K*°) from the event, but significant, and anticorrelated, when selecting single charged 
particles (K^). 
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As a subsidiary check, we replace the atB{PT,nv) parameterization in the likelihood by 
&B(PT, n v) + a ' i where a' is a free parameter. The fit returns a' = —0.001 ± 0.079 for the 
_B°'s, and a' = 0.036 ± 0.052 for the B + result. If a^Qo-r^y) is a good description of the 
data, then a' should be close to zero, as indeed they are. 

The systematic uncertainties of the dilution measurements are determined by successively 
shifting the parameters fixed in the fit up and down by la and repeating the fit. The resulting 
shifts in the fitted dilution are taken as the high and low systematic uncertainties due to 
that input parameter. Section |VI C| discussed the uncertainties assigned to these input 



parameters. The results are shown in Table [XT for J / ipK + and in Table |X1J| for J / ipK 



beginning with the uncertainties associated with the B decay properties. 

Next are the systematic uncertainties arising from the uncertainty in the parameters 
(a's and o's of Eq. (|76|)) describing the signal tagging asymmetry as- Because we used the 
constraint on ■Jb/c^b (App. |C|), we cannot vary as and 7# independently. We vary a B by 
varying the a and b parameters individually by la for fixed "central" values of ^bI&b — 0, 
1, and 2.5, and remaximizing. In this case the dilution shift is from the difference between 
the nominal and shifted values where both use the same fixed value of 'Jb/c^b in the fit. 
The maximum excursion of the dilution among the three 7b /«b combinations is selected for 
each a (o) parameter as the uncertainty for that a (6), irrespective of the ^b/o-b value used 
for the other a's and 6's. The tables list all variations, including those not used. While this 
mixing of ^b/o-b 1 ^ is nominally inconsistent, it provides a conservative estimate. 

The contributions from the uncertainty oyi^b/o^b follow in the tables using the nominal 
as- Table |XII| also includes the effects from the K-tc swap parameters. 

The J I ipK + systematic uncertainty is composed of roughly equal contributions from the 
B lifetime, tagging charge asymmetry, and tagging efficiency asymmetry, but overall has 
a small systematic uncertainty. The largest effects for the J/if)K*° are due to the tagging 
asymmetry and the width of the swapped K-tc mass distribution X$, which has a strong 
asymmetric effect. If the swapped Mm distribution is broad there is little effect; however, as 
it gets narrower it is more difficult to distinguish the swapped from the unswapped events 
and a larger uncertainty ensues. 

The positive (negative) shifts of the dilution due to each parameter are added in quadra- 
ture to obtain the positive (negative) "combined uncertainty" of the dilution in the Tables. 
We thereby obtain the result, 

V+ = 0.185 ±0.052 t fol ( 78 ) 
X?o = 0.165 ±0.112+°;°^. (79) 

These results are similar to the £D^ results of Eqs. (^) and (0). In the next section a 
detailed comparison will be made between the two sets of measurements. 



VII. CHECKS AND COMPARISONS BETWEEN £D^ AND J/^K DATA AND 

MONTE CARLO SIMULATION 

This section presents checks on the robustness of our results and makes a closer compar- 
ison between tagging in our and J/ipK samples as a means of furthering our study 
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of Same Side Tagging. Although the physics processes believed to be responsible for the 
observed flavor-charge correlations should not depend upon the B decay mode, potential 
experimental biases could influence them differently, for instance the fact that the B's in the 
£D^ data tend to have higher p^s than those in the J/ipK data. Given the limited statis- 
tical power of these measurements, we complement the comparisons by also showing some 
results from Monte Carlo simulations. The simulation provides a further means to study 
possible systematic differences between the two data samples, and to gain some insight into 
underlying mechanisms. 

It is not obvious to what extent one can rely upon a given simulation to model particle 
distributions from fragmentation and underlying parton interactions in a pp event. While 
considerable effort has gone into developing and tuning simulations for e + e~ collisions, the 
state of the art is somewhat less well developed for the more complex high-energy hadron- 
hadron collisions. We therefore consider several basic comparisons between data and the 
simulation. We find reasonable agreement and conclude that the simulation is a fair model 
of the data, although our comparisons are necessarily rather coarse. Having developed some 
confidence in the simulations, we proceed to compare £D^ and J/ipK dilutions. It should 
be stressed that prior to this point of comparison the analyses described in this paper have 
not depended upon accurately simulating the tagging (though we have used simulations to 
model decay kinematics, where the models are well established), and indeed, the analyses 
were designed with this independence in mind. 

For the simulations discussed in this section we have used the PYTHIA Monte Carlo 
generator, albeit tuned to match the charged particle distributions in the £D° mode as 
described in Appendix [A 2 b| . Samples of B° and B + mesons were generated for the £D^ and 
exclusive J / ip decay modes. The B decay modes were forced via specific channels for efficient 
generation, and the £D^*' events had a sample composition approximating that found in the 
data. The events were passed through the CDF detector simulation, reconstruction and 
selection code, and finally the SST algorithm. 



A. General comparisons 

We begin our comparison of data and simulation by examining distributions of some 
basic variables. In all these comparisons we use mass sideband-subtracted data samples, 
and we average over the decay modes according to their contribution to the £D^ or J/tpK 
sample. Within each class of sample, the decay-mode specific distributions are very similar 
to one another. 

The multiplicity distribution of SST candidates per B (i.e., tracks which satisfy all 
SST cuts except for the criterion of minimum p™ 1 ) is shown in Fig. |18] for both data and 



simulation. There is general agreement between the simulations and data for both £D^ 
and J/ipK. The £D^*' channel tends to show on average slightly more tagging candidates 
per B than the J/ipK. This difference is due to the higher energies which characterize the 
£D^ candidates and to the extra tracks present from B daughters which are not used in the 
partial B reconstruction. We note as a point of contrast that the corresponding distributions 
for the sidebands show some significant variations across the decay modes, so the agreement 



seen in Fig. 18 is not a trivial result. 
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We impose the full SST criteria and show the pr distributions of the SST tags in Fig. |X9 . 
Again, there is good agreement between the data and simulation for the two types of data. 
The tag px distributions of the two types of modes are also very similar to each other. 
The px distribution for tagging candidates (not shown) have a somewhat harder track p? 
spectrum in the higher-energy ED***' sample than in the J/^K mode, but in both cases the 
data and simulation agree well there, too. 

Finally, we compare the p™ 1 distributions for the tag tracks in Fig. The simulation 
again agrees well with the data, and in this case the £D^ and J/ipK modes are also very 
similar. 

Comparisons using these three variables only provide a limited test, but they are closely 
related to our SST algorithm and indicate that the simulation reproduces basic characteris- 
tics of the data. 



B. Influence of the tagging pt threshold 

Our SST algorithm demands that tag candidates have a minimum pt of 400 MeV/c as a 
compromise between the low-pr tracking asymmetry (Sec. |VI C 3| ) and the declining tagging 



efficiency for an increasing threshold. In this section we consider the influence of this cut, in 
particular, the stability of our results when repeating the analyses for a range of px(SST) 
thresholds. This variable is also a useful vehicle for exploring some of the features of SST. 

Applying a different tagging cut to the same sample means that some S's will be tagged 
by a different particle, others will keep the same tag, and others will no longer be tagged 
at all. If a tag changes, the new tag is largely uncorrelated with the old, so that, with 
respect to the new tags, a statistically independent subset of the data is created. Thus, 
repeating the analysis with different tagging cuts produces measurements which are partially 
correlated with each other. The greater the change in the tagging cut, the weaker the relative 
correlation. We do not attempt to unravel this complex pattern of correlations; we merely 
show variations with the pt cut to show the dependence of the tagging results on changes in 
this parameter. We quote the naive statistical errors from the fits of a pt{SST) scan with 
the understanding that the various points and their errors are correlated in an unspecified 
fashion. 

We first consider the stability of our main physics result, Am^. Figure [H] shows the 
variation of Am^ as the SST p^ cut is varied from 0.3 to 1.6 GeV/c in the ID^*> sample. 
The results are fairly typical of those where the tagging p? threshold is being scanned. 
The values are reasonably stable, i.e., smoothly varying in accordance with the subsample 
correlations mentioned above, and quite consistent with a constant value. 

We also examine the effectiveness of the tagging algorithm as the px threshold is changed. 
Figure ^2] shows the neutral dilutions for the ID^*> and J/ipK*° data. They are both 
relatively constant and are consistent with each other, although the J/ipK*° values offer 
little discrimination. The simulation, also shown in Fig. |22|, agrees well with the data. 



The B + dilutions are shown in Fig. ^3| and display a striking rise as the tagging cut is 
increased. The J /ipK + data, though with sizable statistical uncertainties, show a rise similar 
to £t)W but perhaps offset to lower overall dilution. There is an apparent shape discrepancy 
between the simulation and the J/tpK + data. Appendix [D| describes a x 2 -based comparison 
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between the simulation and data which indicates that statistical fluctuations among these 
correlated measurements can produce such disagreements (or larger) about 22% of the time. 

The different magnitudes and behaviors of the charged and neutral dilutions, as seen 
in Figs. [22| and |23], may at first be surprising in light of the isospin symmetry of the B-ir 
system. However, these differences may result from the fact that tags include not only pions, 
but kaons and protons as well pl[ . For instance, a K~ would be associated with a B + , while 
a B° should be accompanied by a K°, which cannot be a tag. The contrast between the 
charged and neutral dilutions is amplified because when the associated kaon is a K*, the final 
charged kaon is always a K~ , i.e., B°K*° followed by K*° — > K~tc + versus B + K*~ followed 
by K*~ — > K~7i°. A similar argument can be made that the tagging contribution from 
(anti)protons also degrades V and enhances V + . We test this hypothesis in the simulation 
by restricting SST to tag only on prompt pions. The predicted results are shown in Fig. |24] , 
where it is seen that this restriction makes the charged and neutral dilutions nearly identical 
to one another. The ability of the simulation to reproduce the striking behavior of Figs. [22| 
and |23] offers indirect quantitative evidence that tagging on non-pions is the effect in play 
here.[] A similar computation, with the same implications, was reported in Ref. |J for a 
variant of Same Side Tagging in Z Q — > bb; although the measured charged and neutral 
dilutions showed a difference consistent with this effect, the uncertainties were so large that 
no definite conclusion was made there. 

Finally, we show the tagging efficiencies in Fig. |25] as a function of pt{SST). We again 
average all the modes of a given type. The data do not show a clear difference in the 
efficiencies among the separate modes, but the simulation indicates that the efficiency for 
charged B mesons is shifted higher than the neutrals by ~ 2% (absolute) consistently over 
this pt(SST) range. The calculated shapes agree fairly well with the data, but in the case 
of £D(*' at least, there is a small systematic shift in the efficiency: in the simulation there 
are too few cases where there is no SST candidate associated with the £D^*\ This effect 
can also be seen in Fig. where the simulation is slightly below the data in the zero bin. 
However, the shape of the efficiency curve tracks the data well in Fig. [25|, reflecting the good 
description of the track px distribution. Also note that the efficiency falls off much more 
slowly than the tag pt distribution (Fig. [L9|), because as the threshold is raised and tags 
fail the px cut, another track will often be available as a tag. The J/ipK data shows some 
tendency to have a higher efficiency than the calculations, but the difference is statistically 
marginal (recall that the points are correlated and fluctuations manifest themselves over a 
range of bins) . 

As noted in Sec. [TTT| , the error on an asymmetry measurement scales with the "effective 
tagging efficiency" eV 2 [Eq. (0)]. The decline in efficiency with increasing p T (SST), along 
with the relatively constant neutral dilution, imply that the optimum SST threshold should 
be relatively low for B 0, s and somewhat higher for B +, s. Our a priori choice of 400 MeV/c 
was a balance between the falling efficiency as the pt threshold is increased and the inherent 
tracking asymmetry at low p? (Sec. [VI C 3| ) . These simulations indicate that the maximum 
eD\ occurs for a px threshold of about 400-500 MeV/c for our £D^ samples, and about 



9 This observation indicates that the SST dilution for neutral -B's can be significantly improved 
with particle identification. 
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600-700 MeV/c for the J/ipK*°. The rising charged dilutions compensate for the falling 
efficiency such that eT>^_ shows a broad plateau starting at about 600 MeV/c for the £D^*\ 
For the J/i/jK + , eD+ reaches a maximum around 700 MeV/c, but then declines for large 
Pt(SST). Except for eT>^_ from £D^*\ the data are not sufficiently precise to confirm these 
predictions for the optimum threshold. Although our SST threshold was not based on this a 
posteriori analysis, the pt cut is in fact close to the optimum suggested by the Monte Carlo 
simulation. 



C. Dilution comparison between £D^ and J/ipK data 



A cursory comparison of Eqs. ( p7|) and flBq) to Eqs. ([781) and (|79|) already shows that the 
dilutions measured in the £D^ and J/ipK samples are very similar. We consider here how 
well they should agree. 

The main difference between the £D^ and J/ipK data samples lies in their different 
Pt(B) ranges. The £D^ sample requires a single-lepton trigger which has a higher lepton 
Pt threshold than the two-lepton trigger used in the J/i/jK samples. The average pt of the 
B mesons in the £D^*' sample (based on the corrections of Sec. |VC|) is about 21 GeV/c, 
whereas it is about 12 GeV/c in the J/ipK samples. The spectra are shown in Fig. 

We look for a Pt{B) dependence by dividing the data samples into Pt{B) bins and 
repeating the analysis separately for each bin.[^] The results are shown in Fig. No 
apparent Pt(B) dependence is observed, though the statistical sensitivity of the data is very 
limited. The dilution from the J/ipK + point around 15 GeV/c is anomalously low, but the 
other measurements are consistent with the £D^ values, suggesting that the low point is 
simply an unusually large fluctuation. 

Since the data samples are too small to be sensitive to a Pt{B) dependence in the dilution, 
we turn to Monte Carlo simulation. We again use our tuned PYTHIA generator; however, 
in order to generate the very large samples needed for an accurate study we dispense with 
the full detector simulation and instead make simple fiducial cuts and apply a p^-dependent 
track efficiency parameterization. We remove all non-prompt particles from consideration 
as tags in lieu of the SST impact parameter significance cut. The dilutions calculated from 



this simulation are also shown in Fig. ^7], and they exhibit a common shape, rising with 
Pt{B) up to about 15 GeV/c, above which they fall slowly. The ratio of charged to neutral 
dilutions, around 1.35, shows no significant dependence on p T (B).f^ 

We can use the calculated pr(-B)-dependent dilution along with the Pt{B) spectra from 
data to compare the data and simulation without having to subdivide the data into even 
smaller subsamples, as was done for Fig. |27|. The pr(i?)-weighted average dilutions appro- 
priate to each data sample are shown in Table [XIII| . The simulation reproduces the data 



measurements quite well. We also calculate the ratios between the data and simulation 



10 In the case of the £D^ analysis, we now fix Am^ to the world average as in the J/ipK analysis. 

n The dilutions are also observed to be insensitive to the B pseudorapidity, where acceptance and 
trigger effects could be important. 
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values; we find that the ratios are all consistent with one another, and also with 1.0 (the 



measured J/ipK + dilution being 1.3a low; see Fig. |27|) . 



D. Extrapolating dilutions 

The fact that the simulation agrees well with the data, as exemplified by the comparisons 
between the calculated and measured dilutions in Table |Xill| , suggests a method by which 



the SST dilution in any B sample of a similar p? range can be determined. Such knowledge 
is essential for measurements involving B tagging where, unlike in the B° -B° mixing case, 
the dilution is not given by the analysis itself. 

The average dilution for a given sample of B's is calculated by weighting the Monte 
Carlo dilution shape by the sideband-subtracted Pt{B) spectrum for that sample, as above. 
The dilution extrapolation is obtained by multiplying this average by the factor V data I'D mc- 
Since the simulation describes both the neutral and charged dilutions well, one can incorpo- 
rate both of these dilutions in determining the ratio T> data /V MC . This factor is 0.906 ±0.101 
when averaging all the ratios in Table |Xill| (including the correlations between the £D^ 



measurements) . 

The uncertainties on such a dilution extrapolation come from both the measurement 
uncertainties of the £D(*> and J/ipK dilutions (shown above) and from the modeling uncer- 
tainties of the simulation. To estimate the latter, we vary the parameters that control the 
simulation over a fairly wide range. The variations used for the tuned PYTHIA are described 
in App. |A 2 4 We note, however, that the Monte Carlo-derived dilutions always enter the 



above calculation in ratios, i.e., the relative variation of £D^*' to J/ipK, and charged to 
neutral, dilutions as a function of pt{B). We have studied the variations of these ratios 
as we change the inputs to PYTHIA. The largest change in the ratio of £D^ (high-p^) to 
J/ijjK (low-pr) dilutions is 8% from changing the fragmentation p? width to 360 MeV/c. 
The ratio of charged to neutral dilutions shifts by at most 4%, also when the fragmentation 
p T width is set to the low value. 

We have seen (Fig. |4|) that the dilutions are also affected differently by tagging on 



pions and non-pions. This difference introduces another source of uncertainty, especially 
when relying on T> + to constrain T>q as suggested above. One could forgo this additional 
constraint and accept a somewhat larger error on T>o, but the extrapolation using T> + and 
V is not unduly sensitive to the fraction of non-pion tags. As discussed in App. |A 2 cj , we 



estimate the uncertainty due to this effect by allowing the K + to 7r + ratio of tags to vary by 
±30%, and the p to ir + ratio by ±50%, and find that the ratio V /V + changes by ±0.084 due 
to the kaons and ±0.045 for protons. A simple extrapolation from a T> + measurement to T>q 
would translate into a neutral dilution uncertainty of ~ 0.02 (given a V + of 0.25). However, 
since we propose extrapolating from both charged and neutral dilution measurements to a 
neutral dilution, the species-sensitive scaling factor (V /V + ) only applies to the charged 
measurement. The significance of the uncertainty on the fraction of non-pion tags is thereby 
reduced in this application. 

As an example, if we applied the above prescription to calculate the dilution appropriate 
for our J/i/jK*° sample using all four measurements of Table |Xill|, then the dilution obtained 



would be V = 0.171 ± 0.019 ± 0.013 [34]. The first uncertainty is from the statistical 



uncertainty on the scale factor above. The other is the systematic uncertainty from the 
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extrapolation. The latter uncertainty is the quadrature sum of the uncertainties obtained 
from varying the PYTHIA parameters (±0.010) and from varying the fraction of non-pion 
tags (±0.008). The systematic uncertainty from the non-pion contribution is one of the 
larger uncertainties, and could be largely eliminated by not using the V + measurements. 
This, however, would result in an overall larger uncertainty due to the increased statistical 
error. On the other hand, the systematic uncertainty is relatively small even using T> + , and 
the dilution determination is not very sensitive to the simulation. 

This general method can be used to estimate the dilutions in a variety of B meson 
samples of interest for precision measurements of CKM parameters that will be performed 
in the future. In the upcoming Run II, CDF expects to have tens of thousands of exclusively 
reconstructed B° and B + decays through various channels, including B° — > J/ipKg, im, 
Kir, and B°' + — > D^n. The above recipe can be applied to all of them in spite of their 
likely differences in selection criteria. Moreover, further dilution measurements can easily be 
incorporated into T> data /V M c, thereby facilitating increasingly precise measurements. The 
individual large exclusive samples will yield good dilution determinations relatively quickly 
and easily, and refined determinations combining different modes can follow. 

E. Comparison Summary 

We have made several checks on the tagging characteristics of SST in our data. The £D^ 
measurement of Arrid is largely insensitive to the pr threshold of the tagging algorithm. How- 
ever, the charged dilution shows a dependence on the pr threshold of the tagging algorithm, 
in contrast with the largely threshold-independent behavior of the neutral dilution. This 
pattern is apparent in both £D^ and J/tpK analyses, and is also reproduced by the Monte 
Carlo simulation. The simulation indicates that this difference is due to tagging on charged 
kaons and protons. 

Comparisons with other variables show considerable consistency in the characteristics 
of the tagging across decay modes and with Monte Carlo simulation. The general agree- 
ment with our tuned PYTHIA is good, although the current data samples are insufficient 
to confirm some of the more subtle behavior suggested by the simulation. In particular, 
the tagging does not appear to be particularly sensitive to the different kinematics of the 
exclusive and semi-exclusive samples in the variables examined, such as the pr of the B, 
suggesting that SST has more general applicability than to the decay modes examined here. 

VIII. SUMMARY 

We have developed a Same Side Tagging (SST) method based on the flavor-charge corre- 
lations between B mesons and nearby charged particles ("7r ± ") at production. We have used 
SST to tag the initial B flavor in two classes of B reconstruction, while the (nominal) decay 
flavor was obtained from the B reconstruction. Comparison of the initial and decay flavors 
allows one to quantitatively study the strength (i.e., dilution) of the B-ir ± correlation, and 
to observe the B°-B° flavor oscillations. 

The first sample consisted of B — > £D^X decays, which, because it is only a par- 
tial reconstruction, involved additional complications. We have discussed extensively these 



45 



complications, including the separation of B + and B° decays as well as the corrections for 
tagging on B decay products. We observed B-n ± correlations, used them to reveal the 
time-dependent flavor oscillation of B^s, and measured its frequency to be 

Am d = QA71±° -%H ± 0.034 ps" 1 . (80) 

This result is comparable to other single tagging measurements, and agrees well with a 
recent world average of 0.484 ± 0.026 ps _1 from 6 measurements |29| 



The dilution of the flavor-charge correlation for this £D^ sample was found to be 0.27 ± 
0.03 ± 0.02 for the charged, and 0.18 ± 0.03 ± 0.02 for the neutral mesons. The effective 
tagging efficiencies eD 2 are 

eD 2 + = 5.2 ±1.2l°;jj% (81) 
tD 2 = 2.4 ± 0.7±g:i %, (82) 

which are the largest values demonstrated to date for tagging methods applied to high 



energy hadron collider data ||18|| . Although other tagging methods may actually have higher 
dilutions than SST, the combination of good dilution and very high tagging efficiency for 
SST results in the largest eD 2 . 

This SST method was further tested in the exclusively reconstructed B + — > J/i/jK + 
and B° — > J/ipK*° decays. The flavor-charge correlations were observed, and the flavor 
oscillation was again seen with the B 0, s; however, the small sample size did not permit an 
accurate determination of Am^. The dilutions measured in these samples agree well with 
those obtained from the £D^ data, although with much less precision. 

The behavior of SST was also studied by comparing the two classes of data samples 
to a version of the PYTHIA Monte Carlo tuned to charged particle distributions from our 
£D° data. Comparing the behavior of several kinematic quantities, the data and the simu- 
lation both portray a consistent picture, indicating that the simulation captures the basic 
features of this SST. Of particular note, the differences in the charged and neutral dilutions 
are principally due — according to the simulation — to tagging on kaons. Also, despite the 
different kinematics of our £D^ and J/ijjK selections, the tagging largely behaves the same 
way for both, as exemplified by the weak dependence of the dilutions on the pt of the B. 
Furthermore, we have developed a general method to estimate the dilution in a sample of 
B mesons starting from the dilution measurements in £D^ and J/ipK samples. We also 
expect this dilution to be fairly close to the dilution observed in £D^*\ provided that the 
average B momentum is not vastly different. 

This Same Side Tagging method has been demonstrated to be a powerful means to tag 
the initial B flavor, even in the complex environment of a hadron collider. In the upcoming 
Run II of the Tevatron we expect to collect ~ 2/6" 1 of data with the upgraded CDF 
detector. This should result in tens of thousands of exclusively reconstructed B° and B + 
decays in various channels that can be used for precision measurements of CKM parameters. 
The Same Side Tagging technique will be useful for those measurements where initial flavor 
determination is critical. 
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APPENDIX A: MONTE CARLO SIMULATIONS 

In this paper, two types of Monte Carlo simulations are used. Calculations depending 
only on the production and decay of B mesons employ a Monte Carlo generator that sim- 
ulates only a single B. Situations which depend upon the fragmentation particles resulting 
from the hadronization of the b quark, as well as the "underlying event" particles, use the 
full event generator PYTHIA. 

The entire £D^ analysis uses the single B generator simulation, with the one exception 
of the determination of the ^c(^) shape (Sec. |V E), which uses the default PYTHIA 



simulation. The J/ ipK analysis also uses the single B generator, and the comparisons made 
in Sec. |V11| rely on a specially tuned variant of PYTHIA. 



1. Simulation of a single B meson 



Monte Carlo simulation of only a single B meson is based on the following elements. 
Single b quarks are generated using the inclusive 6-quark production calculation of Nason, 
Dawson and Ellis [35|, and the MRSDO |36j parton distribution functions. The b quark 



is then transformed into a B meson, with no additional hadronization products, using the 
Peterson fragmentation model (e = 0.006) |T|. The B meson is decayed using the QQ 
program (Version 9.1) |38| developed by the CLEO Collaboration. The sample composition 



parameters governing the B decay are listed in Table XIV 



2. Monte Carlo simulation of the whole event 



a. "Default" PYTHIA 

The PYTHIA Monte Carlo (PYTHIA 5.7/JETSET 7.4) || is used in instances where 
more than just a single decaying B meson is required. PYTHIA simulates a complete pp 
interaction: the bb pair, the hadronization products, and the remaining beam fragments 
("underlying event"). PYTHIA uses an improved string fragmentation model tuned to 
experimental data, mostly from high energy e + e~ collisions. 

Our PYTHIA generation uses most of the typical default parameters. The CTEQ2L EC 



parton distribution functions are used, and the b quarks are fragmented using the Peterson 
fragmentation model (e = 0.006) ||37|| . B** states are also generated by the fractions listed 
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in Table pCVj However, we suppress the actual B decay performed by PYTHIA and instead 
invoke the QQ program with the same parameters in Table XIV| . In this way we maintain 
a consistent decay model across the two different generators. 



b. "Tuned" PYTHIA 



The PYTHIA generator is controlled by a series of parameters whose default values 
have been adjusted to achieve good agreement with, primarily, high energy e + e~ data. 
Discrepancies between the "default" PYTHIA (as defined above) and CDF pp data are 
apparent, especially when considering particle production that does not originate from the 
b hadronization, i.e., the "underlying event." We have made a separate study [41] of the 
fidelity of the "default" PYTHIA generator (after detector simulation) by comparing it to 
the £D° data (Sec. [V A| ). This comparison studied track multiplicities (with SST quality 
cuts) in AR and A0 intervals around the B direction, and in several pt bins. The data are 
found to have a higher multiplicity of underlying event tracks (as measured away from the 
B, e.g., AR > 0.4) than PYTHIA predicts. 

We may obtain a good description of the charged particle multiplicities and px distri- 
butions by adjusting several PYTHIA parameters. The properties of multiple interactions 
and beam remnants are controlled primarily through the multiple interaction cross section 
[PARP(31)], the model for their generation [MSTP(82)], the ratio of gg and qq multiple 
interactions [PARP (85,86)], and the width of the Gaussian pt spread of particles produced 
in the breakup of color strings [PARJ(21)]. Once these parameters are adjusted to obtain 
agreement with the data away from the 6-jets, we assume the underlying event is well mod- 
eled. We then adjust the Peterson constant PARJ(55) so that the generated multiplicity of 
tracks inside the |Ai2| < 1 cone around the b matches the observed one. Table [XV1| lists the 
default and tuned values of the relevant PYTHIA parameters. More details may be found 
in Ref. g|. 

As an example of the effects of the tuning, we show in Fig. 28 the pt distribution of SST 
candidates (i.e., tracks that satisfy the SST selection cuts except for the p™ 1 requirement 
[Sec. [IV B[ ) in £D^*' data and the two simulations. The tuning procedure uses this distri- 
bution from the £D° subsample — except that the tracks were not restricted to AR < 0.7 
around the B as they are in Fig. |28] — so the agreement of the tuned version with the data is 
to be expected. The shape of the default PYTHIA spectrum shows a clear disagreement 
with the data, with a large excess of tracks at low p?. While there is much better agreement 
between data and the tuned PYTHIA in the shapes of the pt and the frequency distribution 
of SST candidates (see Fig. [TH]) , the tuned Monte Carlo underestimates the number of £D^*> 
events that fail to tag by a few percent (see Fig. |25[). The charged and neutral dilutions as 
a function of the SST px threshold (similar to Figs. 21 and 22) for the two versions of the 
simulation differ by not more than ~ 2a of the Monte Carlo statistical uncertainty; despite 
the better description of the data by the tuned PYTHIA, the tagging results are not much 
different between the two simulations. 

We then have a variant of PYTHIA tuned to our £D° data. This was done for only 
one £D^ mode, and comparisons with the others, or the J/ipK modes, are independent of 
the tuning. We find the tuned version generally provides better agreement than the default 
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version with all the data samples considered in this paper, in spite of the fact that the tuning 
used only global multiplicity and pt distributions and did not consider particle correlations. 



Systematic uncertainties for the dilutions derived from PYTHIA 



For the study presented in Sec. VII D| , we rely on the tuned PYTHIA to calculate the 



dependence of the dilution on the px of the B meson being tagged. In order to determine a 
systematic uncertainty on the dilution extrapolation due to the simulation, we regenerated 
Monte Carlo samples varying selected PYTHIA input parameters. 

The four parameters we varied were the string fragmentation model parameter o^™ 9 
which describes the distribution of particle momenta transverse to the string direction, the 
underlying event cross section scale factor [PARP(31)], the Peterson fragmentation param- 



eter e&, and the combined contribution of the B** modes (see Tables jXV] and |XVj). These 
four were selected as parameters that most directly influence the track momentum and 
multiplicity of potential tags, and hence the dilution. 

We varied vl™ 9 from our tuned value down to 0.36 (the default value) and up to 0.8, 



though the statistical uncertainty from tuning this parameter on the data was only 0.02 fl41 
Likewise, we varied the cross section scale factor from 1.0 (the default value) up to 2.5, 
even though its tuned uncertainty was only 0.04. The large ranges we used for these two 
parameters were chosen as conservative allowances for the applicability of this model.0 The 
range was selected by varying the parameter to (approximately) span a symmetric range 
about the tuned value that included the default PYTHIA value. 

The parameter and the fraction of B mesons originating from B** have been measured 
elsewhere, and their effects on the model are better understood. We varied from 0.004 
to 0.008 and the B** fraction up and down by 25%. These ranges are indicative of the 



statistical uncertainties derived from the tuning studies JIT . 

As an additional systematic uncertainty on the behavior of the dilutions, we varied the 
fractions of kaon and proton tags in our samples. Section |VllIj| indicates that the difference 



between T> + and T> in the simulation is due to tagging on kaons and protons. We varied 
the kaon fraction by ±30% and the proton fraction by ±50% to evaluate this uncertainty. 



3. Detector simulation 

The outputs of the physics simulations are passed through the standard CDF fast detector 
simulation. This simulation is based on parameterizations of detector responses determined 
from data, often test beam measurements. The detector simulation output can be recon- 
structed using standard CDF software. These reconstructed Monte Carlo events may then 
be treated as real data in the analysis. 



12 The PYTHIA default, 

a vr g — 0-36, results from tuning to LEP data. This "string-breaking" 
parameter should be, to first order, the same for e + e~ and pp colliders. The sizeable difference 
with the LEP value may signal a limitation of the tuning procedure, or be a hint that the model 
is inadequate. 
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The inclusive lepton trigger introduces a strong kinematic bias in the £D^*> analysis. 
This bias must be well modeled in the simulation to obtain the proper relative reconstruction 
efficiencies and ct corrections, otherwise an incorrect sample composition will result. We take 
an empirical approach rather than simulate the trigger directly. The trigger is modeled by 
a simple error function parameterization of the ratio of the observed lepton px distribution 
in the data to that generated by the simulation |[27|| . Examples of such ratios and the error 
function fits are shown in Fig. |29] for one signature. Only the region < pt{£) < 20 GeV/c 
is fit, since this is where the effect of the trigger turn-on is the most pronounced. Fits 
are performed on all five decay signatures, and the sample-weighted average of the five 
sets of fit parameters is used to describe the electron and muon trigger efficiencies. These 
parameterizations are then applied to the Monte Carlo events to obtain simulated data sets 
with the correct trigger turn-on. 

A comparison of some kinematic distributions from the data and the simulation is given 
in Fig. ^ for a sample decay signature. As can be seen, this procedure provides a fairly 
accurate representation of the data. 

For the J/i/jK we found that the Pt{B) distribution that results after detector simulation 
and selection cuts compares fairly well with the data without additional trigger simulation. 
The trigger turn-on at very low Pt(a0 is largely governed by the energy loss in the material 
before the muon chambers. This effect is already included in the detector simulation, and 
thus no specific trigger simulation is done for the J/ipK Monte Carlo samples. 



APPENDIX B: CHARGE ASYMMETRY TAG CORRECTIONS FOR THE J/ipK 

SAMPLES 



A charge bias in the SST algorithm could fake an asymmetry. The maximum likelihood 
approach described in Sec. |VI B| provides a natural tool for the parameterization of a charge 



bias and its effect on data. 

A tagging algorithm is characterized by the probability, P(s\p), that a given production 
flavor p yields a tag of charge s. The production flavor p follows the same convention as the 
reconstructed flavor r: p = +1 for B + and B°, and —1 for B~ and B°. The tag, s, takes 
the value of +1 for tagging on a positive track, —1 for a negative track, and if there is no 
tag. This probability can be written in a form similar to the expression of P(s\r) used in 
our likelihood function (e.g., Eqn. (|70|) ), namely 

pMp)=( 1+pk * (3)V * )us), (bi) 



which is characteristic of an asymmetry in pn^s) with amplitude (<f) = S, P, or L for 
the type of event). 

The six P^slp) describing a tagging method (for a given 0) are reduced to four by the 
two constraints 



p^+Ip) + P<t>{-\v) + PMv) = 1 (B2) 

for either p. Four independent variables may be chosen to describe the tagging in terms of 
these probabilities as 



50 



p 


1+ 


1 \ 
-J- j 


_|_ p 


i _ |_ 1 


P f 


\\ 
~T t 






1 + 


+ , 


I P 

+ -r </> 


| 




i 

+, 


_l r> ! _i_ l "\ 


p 


1 + 


T J 


-1- P, 


l + I J 


-I- P, f 


i ^ 














2 








( + 


+ ) 


+ P 


(+I-) 


~^(- 


+) 


-^(-h) 




( + 


+ ) 




(+h) 


+ ^(- 


+) 


+ ^(-h) 




( + 


+ ) 




(-1 + ) 






-P<t>{-\-) 




( + 


+) 








+) 


+ P*(r\-) 



n, = ;:" : : — i ' — ! ' — :r • ; : (B3) 

= ■ - ■ — (B4) 

(B5) 

(B6) 



The first quantity is the usual dilution,0 and the second quantity is the charge- 
averaged tagging efficiency The charge bias in the tagging algorithm is given as [Eq. ([75])], 
and 5$ is the flavor asymmetry in the tagging efficiency [Eq. ([77])]. We find it more convenient 
to express the latter asymmetry as 7^ = 

Solving for P^(s\p) in terms of V^, e^, a^, and 7^, and then casting the expressions in 
the form of Eq. (|B~T|), one derives the following expressions for n<f,(s), the charge asymmetry 
corrected tag, and £</>(s), the corrected efficiency: 



SJs) 



ttS^) for s 

for, 



^(1 + sa^) for s = ±1 
2(1 - e^) for s = 0. 



±1 




(B7) 
(B8) 



Notice that the untagged events may actually have a small finite dilution since ^(0) need 
not be zero. This non-zero dilution arises because the untagged events contain a greater 
number of events which should have been tagged with the sign against which the tagging 
efficiency is biased. 

These equations provide us with a formulation to incorporate tagging asymmetries in 
the likelihood function (Sec. |V1B|) . 



APPENDIX C: CONSTRAINTS ON THE TAGGING EFFICIENCY 

ASYMMETRY 



As discussed in the latter part of Sec. |VI C 3| , a tagging method may not tag on b and b 



mesons with equal efficiencies. The efficiency asymmetry 5^ is given by Eqs. ( [77] ) or ( |B6[ ), 
and appears in the likelihood function via 7^ = S^/V^ in Eq. (|B~T|). We determine Sp t i for 
the J/ipK backgrounds by letting them float in the likelihood fit. However, for the B signal, 
we independently constrain 7^ as explained here. 



13 With this convention, the dilution is positive for tagging B 0, s and negative for B +, s where the 
sign correlation is reversed. However, we explicitly invert the sign in front of the "rKe(,s)2?£" term 
in Eq. (70) for J/i))K + , so that T>b is positive in Table [X| for both B° and B + . All the background 
dilutions follow the nominal convention of Eq. (|B3|). 
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With an ideal detector the tagging method would be described by some "true" dilution 
and efficiency, and the a and S asymmetries would be zero. A detector bias could alter this 
situation by adding or losing tracks based on their charge. For example, positive tracks may 
be added to the event by proton spallation from the beam pipe. This effect is actually very 
small, but in any case, it adds tracks equally around both b and b mesons. This generates a 
non-zero a (more positive tags than negative) in what was an ideal detector, but 5 remains 
zero (b and b mesons have the same positive tag excess). On the other hand, preferential 
loss of one charge makes 5^0. An efficiency asymmetry is created since B° events are more 
likely to tag on negative tracks than B 0, s — this correlation, after all, is why SST works. The 
CTC has such a reduced efficiency for low px negative tracks (Sec. |V1 C 3| ). 

We consider the situation where we have a net loss of negative tracks, as is actually 
observed in our data. Losing a track has one of three outcomes. First, if there was no other 
SST candidate in the event, the tag would simply be lost, giving a net positive tagging 
asymmetry. If, on the other hand, the SST tagged on another negative track, then the loss 
has no effect, since it is only the sign of the tag which matters. However, if the SST tagged 
instead on a positive track, then the tagging asymmetry would be enhanced over that from 
simply losing negative tags. 

If and are the nominal tagging dilution and efficiency in the absence of negative 
track loss (a' = 5' = 0), then the probabilities P(s\p) [Eq. (|Bl| )1 with the negative track loss 
can be rewritten in terms of the nominal quantities as 

(C1) 



p(+\+) 




PH+) 


=< 




=< 


p(-H 


=4 



2 J "V 2 
l-V' 



2 



(C2) 

Vh (C3) 



l + V' 



(l-Ty), (C4) 



where rj is the fraction of negative tags which are lost but not counting those which re-tag 
on another negative track, and /i(2) is the fraction of p = +1 (p = —1) [or B° (B )] which, 
having lost a negative tag, re-tag on a positive track. 

We can calculate the ratio by substituting Eqs. (|C1|) -(|C4D into Eqs. (|B3|), (|B5| ), 



and (|B6|) , and obtain 



{2 -7/(1 - / + P;A/)}{^[! -f]+Af} 
{1 + 7 - V'Af}{2V' - V (V'[1 + J] - A/)} ' 



(C5) 



where f = (fi + /2V2 and Af = (fi — f^)/^- The behavior of ^b/c-b is shown in Fig. \TL 
for values of a# = 2% and V' B = 16.5%, values which are close to what is observed in data. 
We also use 

V (l + ^) + (l-^)(/-D'A/)' 
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obtained from Eq. (P5|). The nominal dilution used is actually the observed dilution in 
J/ipK data after negative track loss (Section |VI Q ), but the two dilutions are expected to 
be similar in light of the small charge asymmetry of ~ 2% (Sec. |V1 C 3| ). 
The largest j^/a^ is achieved with / = Af = 0.5, which gives 

70 = (P; + 1)(4-(1 + P;)7y) 

a<t> (3-Py(4^+(l-3P^)- 1 > 

For values of and D' B close to what is expected in signal events, this maximal value is 
about 2.5. The maximum corresponds to the unrealistic situation where all the lost negative 
tags (which do not re-tag negative) always re-tag on a positive track for B 0, s and never re- 
tag on a positive track for B 0, s. For the likelihood fit we choose the nominal value of this 
ratio to be 1.0, and for the purpose of evaluating systematic uncertainties this ratio is varied 
between and 2.5. Since as is on the order of 2%, j B lies between and 5%. 



APPENDIX D: STATISTICAL SIGNIFICANCE OF T> + VS. P T (SST) SHAPE 

DIFFERENCES 



The variation of the J/ ipK + charged dilution versus the cutoff shown in Fig. |23| appar- 
ently does not agree very well with either the £D^ data or the tuned PYTHIA simulation 
around pr(SST) ~ 0.6 GeV/c. As noted in Sec. |VIIB| , neighboring points are highly cor- 



related and it is difficult to judge the significance of trends across several points from the 
drawn error bars. The correlation is complicated because events which lose their tags as the 
Pt(SST) cutoff is raised will sometimes re-tag on another, higher p^, track in the event. 
This effect causes the fluctuations to be larger than might be naively expected. We consider 
here a test to gauge the statistical significance of the shape differences. 

The dilution differences between adjacent pt(SST) cut-offs (say from 0.4 to 0.5 GeV/c) 
are much less correlated than the absolute dilutions. The common components largely cancel 
in the differences. We estimate the probability to obtain dilution differences similar to the 
J/ipK + data and use that estimate as a measure of the statistical likelihood to obtain shape 
disagreements like the data. 

We calculate a \ 2 comparing the dilution differences between the data and the Monte 
Carlo simulation, i.e., 



c 2 = E 



Si — 5i x 



a(5i 



;di) 



where i is the index of the p^-cut (13 values in 0.1 GeV/c increments starting from 0.3 
GeV/c), 5i = V i+ i — T>i, T>i is the measured dilution at i, 5i is the corresponding difference 
from the tuned PYTHIA (Sec. [A 2 b|) , and a(5i) is the statistical uncertainty on 5i. We 



calculate these differences relative to the PYTHIA value 5i since T> + varies with pt(SST) 
(Fig. |23|) and would otherwise introduce an unwanted systematic contribution to ( 2 . 



We subdivide the PYTHIA sample into 100 subsamples, each with statistics equivalent 
to the J/ipK + sample, and compute ( 2 for each [Q. These subsamples should have the 
same sort of statistical fluctuations as the data. The ( 2 distribution for the Monte Carlo 
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subsamples is shown in Fig. [32[ The value obtained from the J/ipK + data, ( 2 = 16.7, is 



marked by the vertical line. The data is higher than typical, but well within the spread of 
the Monte Carlo samples. 

The distribution of the ( 2, s should, if the dilution differences are truly uncorrelated, follow 
the standard y 2 distribution for n degrees of freedom, which in this case is the number of 
differences. A fit of the ^-distribution to the 100 subsamples, with nasa free parameter, is 
also shown in Fig. [I2|. The fit yields n = 12.71 ± 0.48, in good agreement with there having 
been 13 differences in the ( 2 sum. 

We compute from this fit the probability for a sample the size of the J/if)K + data to 
yield a ( 2 at, or above, the 16.7 observed in the data to be 22%. Thus, we conclude that 
the observed differences in the dilution shape with pt{SST) are not statistically unusual. 



54 



REFERENCES 



[1] N. Cabibbo, Phys. Rev. Lett. 10, 531 (1963); M. Kobayashi and K. Maskawa, Prog. 
Theor. Phys. 49, 652 (1973). 

[2] UA1 Collaboration, C. Albajar et al, Phys. Lett. B 186, 247 (1987). 

[3] ARGUS Collaboration, H. Albrecht et al, Phys_Lett. B 192, 245 (1987). 

[4] For a recent review of experimental results on B°B° oscillations see O. Schneider, "Heavy 
Quark Spectroscopy, Oscillations, and Lifetime," presented at the 18th International 
Symposium on Lepton-Photon Interactions, Hamburg, Germany, to appear in the pro- 
ceedings. 

[5] HERA-B Collaboration, DESY-PRC 95/01, 1995; CLEO III Collaboration, Cornell 
CLNS 94/1277, 1994; BELLE Collaboration, KEK Report 3-1995, 1995; BaBar Col- 
laboration, SLAC-R-95-457, 1995; DO Collaboration, FERMILAB-Pub-96/357-E, 1996; 
CDFII Collaboration, FERMILAB-Pub-96/390-E, 1996. 

[6] CDF Collaboration, F. Abe et al, Phys. Rev. Lett. 80, 2057 (1998). 

[7] M. Gronau, A. Nippe, and J. Rosner, Phys. Rev. D 47, 1988 (1993); M. Gronau and 
J. L. Rosner, Phys. Rev. D 49, 254 (1994). 

[8] OPAL Collaboration, R. Akers et al, Z. Phys. C 66, 19 (1995). 

[9] ALEPH Collaboration, R. Barate et al, Phys. Lett. B 425, 215 (1998). 
[10] OPAL Collaboration, K. Ackerstaff et. al, CERN preprint CERN-EP/98-001. 
[11] CDF Collaboration, F. Abe et al, FERMILAB-Pub-98/189-E, submitted to Phys. Rev. 
Lett.. 

[12] CDF Collaboration, F. Abe et al, Nucl. Instrum. Methods A 271, 387 (1988). 
[13] D. Amidei et al, Nucl. Instrum. Methods A 350, 73 (1994); P. Azzi et al, ibid. 360, 
137 (1995). 

[14] CDF Collaboration, F. Abe et al, Phys. Rev. D 50, 2966 (1994). 

[15] CDF Collaboration, F. Abe et al, Phys. Rev. D 57, 5382 (1998). 

[16] ALEPH Collaboration, R. Buskulic et al, Phys. Lett. B 284, 177 (1992). 

[17] ALEPH Collaboration, R. Buskulic et al, Z. Phys. C 75, 397 (1997); OPAL Collab- 
oration, K. Ackerstaff et al, ibid. 76, 401 (1997); ibid. 76, 417 (1997); DELPHI Col- 
laboration, P. Abreu et al, ibid. 76, 579 (1997); L3 Collaboration, M. Acciarri, et al, 
CERN-EP/98-028, submitted to E. Phys. J. C. 

[18] CDF Collaboration, F. Abe et al, to be submitted to Phys. Rev. D. 

[19] OPAL Collaboration, R. Akers et al, Phys. Lett. B 327, 411 (1994). 

[20] R.D. Field and R.P Feynman, Nucl. Phys. B 136, 1 (1978). 

[21] I. Dunietz and J.L. Rosner, Phys. Rev D 51, 2471 (1995). 

[22] Particle Data Group, R.M. Barnett et al, Phys. Rev. D 54, 1 (1996). 

[23] For a recent survey of D** production measurements, see C. J. Kreuter, "Measurements 
of B — > D**lv, n presented at the 2nd International Conference on B Physics and CP 
Violation, Honolulu, Hawaii, March 24-27, 1997, to appear in the proceedings. 

[24] D. Ebert, V.O. Galkin, and R.N. Faustov, Phys. Rev. D 57, 5663 (1998). 

[25] CLEO Collaboration, R. Fulton et al, Phys. Rev. D 43, 651 (1991). 

[26] N. Isgur, D. Scora, B. Grinstein, and M. Wise, Phys. Rev D 39, 799 (1989). 

[27] P. Maksimovic, Ph.D. dissertation, Massachusetts Institute of Technology, 1998 (un- 
published) . 

[28] CDF Collaboration, F. Abe et al, Phys. Rev. Lett. 76, 4462 (1996). 



55 



[29] R.M. Barnett et al, Phys. Rev. D 54, 1 (1996), and 1997 off-year partial update for the 
1998 edition available on the PDG WWW pages (URL: |http://pdg.lbl.gov/| ). 

[30] CDF Collaboration, F. Abe et al, Phys. Rev. D 53, 3496 (1996). 

[31] CDF Collaboration, F. Abe et al, Phys. Rev. Lett. 72, 3456 (1994). 

[32] CDF Collaboration, F. Abe et al, Phys. Rev. Lett. 76, 2015 (1996); Phys. Rev. D 54, 
6596 (1996). 

[33] CDF Collaboration, F. Abe et al, Phys. Rev. Lett. 75, 1451 (1995). 
[34] K. Kelley, Ph.D. dissertation, Massachusetts Institute of Technology, 1998 (unpub- 
lished) . 

[35] P. Nason, S. Dawson, and R. K. Ellis, Nucl. Phys. B 327, 49 (1988). 

[36] A.D. Martin, W.J. Stirling, and PC Roberts, Phys. Rev. D 47, 867 (1993). 

[37] C. Peterson et al, Phys. Rev. D 27, 105 (1985). 

[38] P. Avery, K. Read, and G. Trahern, Cornell Internal Note CSN-212, 1985 (unpublished). 
[39] H.-U. Bengtsson and T. Sjostrand, Computer Physics Commun. 46, 43 (1987); 

"PYTHIA 5.7 and JETSET 7.4: Physics and Manual," by T. Sjostrand (Lund U.) 

LU-TP-95-20, Aug 1995. 
[40] H.L. Lai et al, Phys. Rev. D 51, 4763 (1995). 

[41] D. Vucinic, Ph.D. dissertation, Massachusetts Institute of Technology, 1998 (unpub- 
lished) . 



56 



TABLES 



Decay Signatures 



Selection 


Cuts 




£D- 


K+7T- 


ID*~ 

K + ir~TT + TT~ 


K+TT-TT 


pt(£) > 


(GeV/c) 


6.0 


6.0 


6.0 


6.0 


6.0 


Pt(K) > 


(GeV/c) 


0.7 


0.6 






1.0 


p T (n) > 


(GeV/c) 


0.5 








0.8 


Pt(D) > 


(GeV/c) 


2.0 


3.0 








d /a > 




3.0 


2.0 


1.0 


0.5 


1.0 


L xy (D)/a Lxy > 




3.0 


5.0 


1.0 


1.0 


1.0 


\Am(D*)\ < 


(MeV/c 2 ) 






3.0 


2.0 




m(£D) < 


(GeV/c 2 ) 


5.0 


5.0 








-0.5 < ct D < 


(mm) 


1.0 


2.0 


1.0 


1.0 


1.5 



TABLE I. Kinematic and geometric selection cuts for the five decay signatures. The impact 
parameter significance cut do/ao is applied to D daughter tracks. L xy (D)/aL xy is the D decay 
length significance relative to the primary, while ct£> is the proper decay length of the D with 
respect to the B vertex, and Am(D*) is the mass difference between the D* candidate and the D 
candidate plus pion mass. 



Name 


J F 


Width 


Decay Modes 


n* 


0+ 


wide 


Dtt 


D\ 


1+ 


wide 


D*TT 


L>i(2420) 


1+ 


narrow 


D*7t 


£>5(24fiO) 


2+ 


narrow 


Dtt, D*ir 



TABLE II. The expected D** states and properties. The classification into wide and narrow 
states follows the predictions of Heavy Quark Effective Theory. 
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Decay Signatures 



B° Decay Chains £D*~ ID~ £D° 

D**~£ + v 



£)*— _> D*\ 



P* - P°tt° /**iV(2/3)B*(^-7rS) 
P*° -► P° 7 f**P v (2/3) B*(S°7) 

£>- - S°7r- /**Py(l/3)£*(P°7r-)e(7r*) /**iV(l/3)B*CD°7r-)(l - e(ir,)) 

P*~ - £>-7r° f**P v (l/3)B*(D-iT°) 



D*~ -► P~ 7 /**P y (l/3)H*( J D-7) 

£>*-- EV- /**(l-Py)(2/3) 

£>**- _> g-Trg, — /**(!- iV)(l/3) - 

> 

£)*- _> 5V /*fi*p 7r-)e(vr,) /*6*(E%-)(1 - e(^)) 

P*~ -► P^° f*B*{D~^) 

D*~ -► D-7 — f*B*(D~j) — 

I D-t+v — / — 
P + Decay Chains 

£>- - S°tt- /**Py(2/3)£*(P°7r-)e(7r*) f**P v (2/3)B*(D%-)(l - e(vr*)) 

P*~ -► r»-7r° /**Py(2/3)fi*( J D-vr°) 

P*-^P° 7 /**P y (2/3)6*(D-7) 

P* ^!* ^ /**iV(l/3)B*(B°7rO) 

P*° P° 7 /**P y (1/3)23* (P° 7 ) 

S** ^D-7T+ /**(l-Py)(l/3) 

P**°^pV°* /**(!- Py)(l/3) 











P*° - pV° 








P* - P° 7 






/*£*(P°7) 


-► 5Vi/ 






/ 



TABLE III. Table of the various P decay chains and their contributions to the sample compo- 
sition of the three general categories of decay signatures (£D*~ , £D~ , and £D ). Dashes indicate 
that a particular decay chain makes no contribution. The total contribution to a given sample is 
simply the sum over the entries in a vertical column. The branching ratio of P* — ► XY is denoted 
by U B*(XY)." All t^'s, 7r**'s, and 7's are lost from the reconstruction, and 7r^'s are identified 
with efficiency e(7r*). See the text for a discussion of the other parameters. 



58 



Signature 






a ct (0) (um) 

\ J Vr 111 / 




b 


£ + D~, D~ - 


* K+TT IT 




39 




0.108 


l+D*-, D° - 


-» K+TT- 
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0.075 


l+D*-, D° - 






49 




0.073 


l+D*-, D° - 






62 




0.070 


l+D\ D° -> 






45 






TABLE IV. The paramefr 


ers of the linear mo 


del of the ct resolution [Eq. (|53| 


)] for the five direct 


decay chains. 












Decay signature 








J^** (meets) 


l+D~, D~ - 


» i^ + 7T 7T 








0.056 ± 0.022 


l+D*-, D° - 










0.003 ± 0.029 


l+D*-, D° - 










-0.016 ± 0.026 


l+D*-, D° - 


-» ET+vr-vr 








0.034 ±0.021 












0.029 ±0.018 


TABLE V. The fractions of tags (with no (Iq/cfq cut) identified as ir+^ 


candidates, R** , measured 


in the five decay signatures. 










Input 




Fit Input 






Fit Output 


Parameter 


Value 


Error 




Value 


Error 




0.360 


±0.120 




0.309 


±0.100 




2.50 


±0.60 




2.51 


±0.60 


cr (/an) 


468 


±18 




468 


±18 


r+/r 


1.020 


±0.050 




1.021 


±0.049 



TABLE VI. The fit input and output values of the measured sample composition parameters. 
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name 



Am d 


1.000 


-0.372 


-0.172 


v + 


-0.372 


1.000 


0.372 


v Q 


-0.172 


0.372 


1.000 


Rf 


-0.020 


0.126 


0.007 


j ** 


-0.385 


0.406 


0.504 


Pv 


-0.326 


-0.284 


-0.310 


e(vr*) 


-0.031 


-0.082 


0.100 


£,norm 


0.304 


-0.355 


-0.445 


TO 


-0.005 


0.002 


-0.001 


t+/t 


-0.051 


-0.157 


0.009 



TABLE VII. The correlation coefficients for Am^, V + , and Vq, with respect to all ten fit 
parameters. 



Source 


f(Am d ) 


(7(2?+) 


a(V ) 


Sample Composition 
6-quark spectrum 
e isolation cuts 
Decay model 
L xy resolution 
£,Mc{ct) shape 
B -► D { s ] D^X 
B s -► vtDf 
g^cc^ £DW 


+0.0295 
-0.0310 

±0.0060 
±0.0045 
±0.0115 
±0.0033 
±0.0035 
±0.0010 
±0.0010 
±0.0006 


+0.0216 
-0.0131 

±0.0052 
±0.0036 
±0.0005 
±0.0003 
±0.0002 
±0.0006 
±0.0019 
±0.0012 


+0.0225 
-0.0131 

±0.0017 
±0.0047 
±0.0045 
±0.0000 
±0.0015 
±0.0004 
±0.0008 
±0.0025 


Total 


+0.0329 

-0.0343 


+0.0226 

-0.0147 


+0.0237 

-0.0150 


TABLE VIII. Table of the systematic uncertainties. 


Decay 
Signature 


B -> D^D^ 


Fractional Contribution 
B s z^£>** 


g -► £DW 


1+D-, D~ -► K + 7T-7T~ 

l+D*-, D° -► i^T+vr- 

^° -► i^+Tr-vr 


0.017 
0.005 
0.005 
0.006 
0.009 


0.011 
0.008 
0.008 
0.009 
0.011 


0.005 
0.002 
0.005 
0.019 
0.005 



TABLE IX. Fractional contribution of B -► DT'D^X, B -► D^DW, and 5 cc -► £DW 
to the £DW samples. 
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Output Parameters 



B + -► J/ipK + 



B 



o 



J/ipK 



*o 



Frac. Signal 
Frac. L-Lived Back. 
Mass Error Scale 
cr Error Scale 
Prompt: 
Mass Slope (xl(T 3 ) 
Back. Tag Eff. 
Tag Asym. 
Eff. Asym. 
Recon. Asym. 
Dilution 

Long-Lived: 
1st Lifetime (/jm) 
2nd Lifetime (/xm) 
Frac. Neg. Back. 
Frac. 2nd Lifetime 
Mass Slope (xl(T 3 ) 
Back. Tag Eff. 
Tag Asym. 
Eff. Asym. 
Recon. Asym. 
Dilution 

B Signal: 
Tag m. 
Recon. Asym. 
Dilution 



f B 
h 
X 
Y 

Cp 
ep 
ap 
Sp 
Rp 
v P 

TLX 
TL2 

In 
Cl 

£L 
OIL 
S L 

Rl 
V L 

Rb 



0.067 ± 0.003 
0.160 ±0.008 

I. 34 ± 0.05 
0.99 ±0.01 

II. ±0.9 
0.703 ± 0.005 
0.033 ± 0.013 

-0.002 ±0.007 
0.003 ±0.011 
-0.069 ±0.013 

595. ± 53. 
135. ±9. 
0.137 ±0.014 
0.781 ± 0.030 
-12.3 ±2.7 
0.771 ±0.014 
0.015 ±0.037 
-0.026 ±0.018 
0.030 ± 0.034 
-0.089 ± 0.038 

0.624 ± 0.020 
0.077 ± 0.041 
0.185 ± 0.052 



0.156 ±0.009 
0.222 ±0.023 
1.54 ±0.10 
1.06 ±0.03 

4.4 ±2.4 
0.830 ±0.011 
0.092 ±0.031 
0.012 ±0.013 
0.036 ± 0.029 
-0.003 ± 0.031 

371. ±65. 
99. ±21. 
0.096 ± 0.029 
0.626 ±0.104 
-17.6 ±5.3 
0.778 ±0.031 
-0.044 ± 0.079 
-0.029 ± 0.038 
0.095 ± 0.070 
-0.050 ± 0.079 

0.635 ± 0.030 
-0.086 ± 0.068 
0.165 ± 0.112 



TABLE X. Likelihood fit results for B + -> J/t/jK + and B° -> J/i/jK*°. The parameter of 
central interest is the signal dilution in the bottom line. 
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Parameter 


Central 
Value 


Variation 


snirt 

Negative 


in 

Positive 


T~B 




(fim) 


450 


1 1 Q 

±lo 


— U.UUlz 


n nm 1 
U.UU11 






IbIolb = 1 








ai 


x 1U [u-ev /c) 


o.y 


±1.0 


n nnm 
— U.UUU1 


n nnm 
U.UUU1 


I, 


xlO (Gel//cj 


1.3 


1 f~\ A 

±0.4 


n nnno 

U.UUU3 


n nnno 

— U.UUU3 


0,2 


xlO -2 


1.4 


±0.4 


0.0001 


—0.0002 


7 


xl0~ 2 


2.6 


±0.8 


—0.0005 


0.0005 






JB/C^B — U 








a i 


xlU (GeV/c) 


o.y 


±1.8 


U.UUUz 


n nnno 

— U.UUUo 


L 

h 


xlU [Gev/c) 


1 o 

1.0 


i r\ a 

±0.4 


n nnn'7 

0.0007 


—0.0007 




xlU 


1 A 

1.4 


±0.4 


—0.0004 


0.0003 


j 

02 


xlU 


z.b 


±0.8 


—0.0010 


0.0010 






IB 1 OLE = 2.5 








a x 


xlU {LreV/c) 


3.9 


±1.8 


n nm n 
— U.UU1U 


n nnnn 

u.uuuy 


bi 


xKT 3 (GeV/c) 4 


1.3 


±0.4 


-U.U005 


U.U004 


0,2 


x 10~ 2 


1.4 


±0.4 


0.0015 


-0.0019 




xlU" 2 


2.6 


±0.8 


0.0005 


-0.0008 






«B = Central Value 






IB / OL B 




1 


1.5 
1.0 


0.0015 


-0.0034 



Combined Uncertainty 



TABLE XI. The fixed inputs for the B + — ► J/ipK + fit, their central values, la variations, and 
the resulting shifts of the central value of V + . The shifts are combined in quadrature (see text) to 
obtain the combined systematic uncertainty. 
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Parameter 


Central 
Value 


Variation 


onirc 
Negative 


m Uq 

Positive 


TB 




468 


±18 


a a aao 

—0.0002 


a a aao 

0.0002 


Am 


{ps- 1 ) 


0.474 

1b/olb = 1 


±0.031 


-0.0005 


0.0003 


a i 


xlO * (Gel//cJ 


o.9 


±1.8 


A AAA o 

—0.0008 


a AAA 'I 

0.0007 


7 

h 


xlO (Gev//cJ 


1.0 


1 A A 

±0.4 


a nno a 

0.0024 


A AAA V 

—0.0025 


do 


x 10~ 2 


1 A 

1.4 


±0.4 


0.0013 


—0.0014 


b-2 


xl(T 2 


2.6 

7_B / OLB = 


±0.8 


-0.0039 


0.0038 


at 


xlU (Gei//c) 


o a 

3.9 


±1.8 


a nnn o 

—0.0008 


A AAAO 

0.0008 


I 


win— 3 / s~i t /■ / \4 

XlU ° (Gevyc) 


1.0 


1 A A 

±0.4 


a aaa a 

—0.0004 


A AAA A 

0.0004 


a 2 


xlO 


1 A 

1.4 


1 A A 

±0.4 


a nnno 

—0.0002 


A AAA 1 

0.0001 


7 


xl0~ 2 


2.6 

IB 1 olb = 2.5 


±0.8 


r\ AAAO 

—0.0002 


0.0001 


a i 


xlU {Gev/c) 


3.9 


±1.8 


a Am 1 

—0.0011 


A A A 1 A 

0.0010 


h 


xlO (GeV/c) 4 


1.3 


±0.4 


0.0063 


-0.0065 


a 2 


xlO z 


1.4 


±0.4 


0.0040 


r\ r\ r\ a a 

—0.0044 


i. 

02 


X 1U 


2.6 


±0.8 


n nnno 
— U.UUyz 


U.UU88 






as = Central Value 






IB 1 OL B 




1 


+1.5 
-1.0 


0.0031 


-0.0057 


fs 




0.1 


±0.1 


0.0086 


-0.0160 


X s 




5.0 


±2.0 


0.0111 


-0.0003 


MS 




-0.5 


±0.5 


-0.0008 


0.0007 



Combined Uncertainty 



TABLE XII. The systematic uncertainties from the fixed B° — > J/ipK* fit parameters. The 
table is similar to that for J/^K + (Table |Xl|) except for the addition of the oscillation frequency 
Arrid, and fs, an d Xs, which model the K-tt swapping in the K*° reconstruction. The shifts 
in Vq are combined in quadrature (see text) to obtain the combined systematic uncertainty. 



Sample MC Calculation Data Meas. Data/MC Ratio 

B° WW 0.196 0.181 ±0.035 0.923 ±0.179 

B°^J/i;K* 0.189 0.165 ±0.112 0.873 ± 0.593 

B + ^£D^ 0.266 0.267 ±0.037 1.004 ±0.139 

B+^J/ipK + 0.254 0.185 ±0.052 0.728 ± 0.205 



TABLE XIII. Calculated and measured dilutions for the ID^ and J/ipK samples. The 
calculated values are from the simplified simulation. The ratios are of the measured values divided 
by the Monte Carlo calculations. 
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Parameter 



Value 



R f 2.722 

fZ 0.231 

fTon 0.125 

/** 0.356 

Py 0.687 

t+/t 1.014 



TABLE XIV. The values of the sample composition parameters used in the QQ (V9.1) B 
decay program. The resonant (B — > u£D**) and non-resonant (B — > uID^tt^X) fractions, 
and /**„ respectively, sum by definition to /**. 



Parameter Value 

a(B*)/a(B + B* + B**) 0.7625 

a(B** : l Pi)/a{B + B** : 5 = 0) 0.320 

a(B** : 3 P )/a(B* + B** : S = 1) 0.033 

: 3 Pi)/a(B* + 5** : 5=1) 0.099 

: 3 P 2 )/a(B* +B** : 5=1) 0.165 



TABLE XV. The cr(S) ratios represent the relative production rates used in PYTHIA for the 
different B mesons. The relative ratios are labeled by spectroscopic notation or their spin 5, e.g., 
u a(B + B** : 5 = 0)" represents the sum of cross sections for B and B** states with zero spin. 



Parameter 


Default 


Tuned 


Description 


MSTP(82) 


1 


3 


model of multiple interactions 


PARP(85) 


0.33 


1.0 


fraction of color-connected 








gg multiple interactions 


PARP(86) 


0.66 


1.0 


total fraction of gg 








multiple interactions 


MSTP(33) 


No 


Yes 


multiply cross-sections by PARP(31) 


PARP(31) 


1.00 


1.66 


increase cross-sections by 66% 


PARJ(21) 


0.36 


0.6 


°Vt 


MSTJ(ll) 


4 


3 


use Peterson frag, for b, c 


PARJ(55) 


-0.006 


-0.0063 





TABLE XVI. The PYTHIA parameters modified from their defaults in order to agree with 
B i/£D°, D° K+7T- data. 
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FIGURES 
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FIG. 1. A simplified picture of fragmentation paths for a 6 quark. 
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FIG. 2. Schematic drawing of the momentum vectors determining pTp 1 and p r £ l of an SST 
candidate. 
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lepton 




71 

FIG. 3. A typical B — ► z/£D event topology, where a -B meson is produced at the event primary 
vertex, and decays after traveling a short distance into a lepton, a D meson, and a neutrino 
(undetected, and not shown). The D later decays into a kaon and one or more pions. 
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783 ± 21 events 



lep + D* 




Knn Mass (GeV/c 2 ) 
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FIG. 4. The mass distributions of the fully reconstructed D candidates (solid histogram) 
for: £+D*-, D° -► K+tt- (upper left); t+D*-, D° -► K+it-it + it- (upper right); ^+L>-, 
D~ — > i^ + 7T _ 7r _ (lower left) and D° — ► i^ + 7r _ (lower right). The dashed histograms are the 

distributions of the wrong-sign (l^K^) candidates. The number of signal events from the fit (solid 
curve) is indicated in each plot. 
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2515 ± 70 events 




0.135 0.155 0.175 0.195 

Mass(K7m.) - Mass(Kn) (GeV/c 2 ) 



FIG. 5. The distribution of m(Kinr*) — m(KTr) for the signature D*~ — > D°tt~ , D° — > K + 7r _ 7r° 
(7t° not reconstructed). The upper histogram is the distribution of the right-sign ID candidates, 
and the lower histogram is for the wrong-sign combinations. The solid curve is the fit of the 
right-sign distribution, with the fitted background component shown by the upper dashed curve. 
This background shape is obtained from a fit (lower dashed curve) of the wrong-sign data. 




FIG. 6. The state diagram for all possible B — ► ID(*)X transitions. The /, /*, and /** are 
semileptonic B branching ratios to charm mesons. 
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FIG. 7. The distribution of K = p T {£D)/p T (B) (left), and K vs. m{lD) (right) with a fit of 
a quadratic function, for the direct decay £ + D*~, D° — ► i^ + 7r _ . The "error" bars in the /C vs. 
m(£D) plot represent the RMS spread of the /C-distribution in each bin, and not the error on the 
bin mean of K. 



1000 



o 



500 



74% of 126/i 
26% of 633/i 




-0.20 0.00 

ALxy (cm) 



0.20 




A(i//? T r)/0/M 




0.00 

Act (cm) 



o.ooo 
0.10 o.oo 



0.25 



ct (cm) 



FIG. 8. The simulation for the direct chain £D*~ , with D° — > K+tt^: AL xy (B) distribution, fit- 
ted with two Gaussians with relative fractions and RMS values listed (top left); A(1//3t7)/(1//9t7) 
distribution (top right); Act distribution (bottom left); and the RMS of the Act distribution as a 
function of ct along with the linear fit for o-^ d (ctkd) (bottom right). 
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lepton 




FIG. 9. A schematic representation of a B — > vlD** decay. SST candidate tracks originate from 
the primary vertex, while the 7r^ track originates from the B decay vertex. The impact parameter 
of a track with respect to the primary vertex is do while the impact parameter with respect to the 
B vertex is given by cIb- When the B vertex and the primary vertex are well separated, the 7r^ 
track usually has a large do and a small ds- The converse is true for primary tracks. 
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FIG. 10. Monte Carlo calculation of £mc as a function of corrected proper time, ct° d , for decay 
signature £ + D~: no d /<jQ cut (top), and d^/ao < 3.0 (bottom). 
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FIG. 11. Impact parameter distribution with respect to the B vertex (ds/crs) for the £ + D~, 
D~ — > i^ + 7r~7r~ decay signature from the data (no do/ao cut, D mass sideband subtracted). 
Right-sign tags are shown by the solid points. Wrong-sign tags (open circles) are renormalized to 
model the right-sign continuum at large impact parameter significance, as described in the text. 
The right-sign excess near zero is due to the 7r^ tags. 
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FIG. 12. Measured asymmetries as a function of the corrected proper decay length ctkd for the 
decay signatures £ + D° (dominated by B + ), £ + D~ , and the sum of all three £ + D*~ (dominated by 
B°). The three £ + D*~ signatures are combined only for display purposes. The dashed line is the 
result of the fit. 
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FIG. 13. The breakdown of the measured asymmetries of £ + D* into the three £ + D* decay 
signatures "£ + D*-, D° -► K+tt-", "£+D*-, D° -► K+vr-vr+vr-" and D~ -► K + vr-7r~". 

The result of the fit is overlaid. 
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Normalized mass distribution for J/tpK + and J/t/jK*° candidates with reconstructed 
smooth curve is a fit from the full likelihood function of a Gaussian signal plus linear 
parameterization. 
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FIG. 15. Charge asymmetry (a) dependence on the track's 1/pr (top) and number of primary 
vertices ny in the event (bottom) for: Same Side Tags (squares), SST candidate tracks (open 
circles), SST candidate tracks with 6- vertex veto (solid circles), and tracks in a cone away from the 
J ftp direction (triangles). The solid curves are the results of least-squares fits of the asymmetry 
parameterization (see text) to the SST candidates with the 6-vertex veto. 
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FIG. 16. The J/ipK*° lifetime distributions for the signal (top) and sideband (bottom) regions. 
Superimposed on the data are the likelihood fit results (solid line). In the signal region, the B 
component is shown by the dashed line and the prompt plus longdived background by the dotted 
line. 
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FIG. 17. The mass sideband subtracted flavor-charge asymmetry as a function of the recon- 
structed ct (points): left is B+ -> J/ipK+, and B° -> J/tpK*° is on the right. Superimposed on 
the data are the likelihood fit results (solid lines). The insets are scans through the log-likelihood 
functions as the dilutions are varied about the fit maxima. Also shown in the main plots are the 
results of simple \ 2 fits to the points (dashed lines, partially obscured by the solid lines). 




3 4 5 6 7 

Number of SST Candidates 



FIG. 18. Number distributions of candidate tagging tracks in both £D(*) and J/ipK data 
(points) and simulation (shaded bands). The widths of the shaded bands are the statistical errors 
from the Monte Carlo sample size. 
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FIG. 19. The px distributions of the tag tracks in both £D^ and J/ipK data (points) and 
simulation (shaded bands, width indicating the statistical error). 
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FIG. 20. The p™ 1 distributions of tag tracks in both £D^ and J/ipK data (points) and simu- 
lation (shaded bands, width indicating the statistical error). 
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FIG. 21. The extracted value of Am^ from the tagged £D^ as a function of the tag pr 
threshold. Error bars are the naive statistical errors returned from the fit, and they are correlated 
with each other, as are the points themselves (see text). 




0.6 0.8 1 1.2 1.4 1.6 

p T (SST) Threshold (GeV/c) 
FIG. 22. Tagging dilution as a function of tag px cut for B°. Data are plotted with solid circles 
(£D^) and triangles (J/ipK), and the corresponding simulations are shown by the shaded bands 
(width indicates the statistical error). The various points and their errors are correlated with each 
other (see text). 
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FIG. 23. Tagging dilution as a function of tag px cut for i? + data in the ID^ (solid circles) 
and J/ipK (triangles) mode. The corresponding simulations are shown by the shaded bands (width 
indicates the statistical error). The various points and their errors are correlated with each other 
(see text). 
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FIG. 24. Dilution as a function of the tag px cut from Monte Carlo simulation for the J/ipK 
modes of B 0, s (solid square) and B +, s (solid circle). When the tagging is restricted to prompt 
pions only, the neutral dilution becomes the open squares, and the charged the open circles. 
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FIG. 25. Tagging efficiency as a function of the tag pt threshold for £D(*) (solid circles) and 
J/ipK (open triangles) data. The results from the simulation are shown by the shaded bands 
(width indicates the statistical uncertainty). 
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FIG. 26. The sideband-subtracted pr(B) distributions for the (corrected) £D^ (top) and 
J/ipK (bottom) samples. The charged and neutral ID^ spectra are very similar, but the two 
J/vpK spectra are slightly different due to the different cuts on the kaon momenta. 
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FIG. 27. Tagging dilution as a function of B meson px for charged (left) and neutral (right) 
mesons. The dilution measurements are plotted at the pr(-B)-weighted centroid of each bin, and 
the horizontal error bars span the width of the bin (arrows indicate that a bin is unbounded) . The 
closed circles are the £D^ data, and the open triangles are J/ipK's. The dashed curves are Monte 
Carlo calculations of the px dependence, and the dotted lines mark the average dilution for the 
data points. 
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FIG. 28. The px distribution of tracks satisfying SST selection cuts except for the p r ^ 1 require- 
ment: average for iD^> data (solid circles), default PYTHIA (open squares), and tuned PYTHIA 
(open circles). 
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FIG. 29. The ratios of the pt{£) distributions of the data to the Monte Carlo simulation, for 
electrons (top) and muons (bottom), for the £ + D~, D~ — ► K + ir~ir~ signature. The distributions 
are fit with the error function, where pj, is the midpoint pr{£) and a PT the "width" of the turn-on. 
The overall normalizations are immaterial since the £D^ analysis requires only relative efficiencies. 
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FIG. 30. A comparison between the data and the single-S Monte Carlo simulation for the 
decay signature £ + D~ , D~ — ► K + 7r _ 7r _ . The distributions compared are: of the lepton, e and 
/i combined (top left), mass of the ID~ system (top right), the px of the £D~ system (bottom 
left), and pr(£D~) after correcting for the missing neutrino (bottom right). Only the agreement 
in the corrected pr(£D~) distribution relates directly to the analysis. 
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D = 16.6% a = 2% 




FIG. 31. Variation of j^/a^ versus f\ and fa for a<f, = 2% and nominal dilution T>', = 16.6%. 
Note that (0,0) is the leftmost corner and (1,1) the rightmost, and that the function is not defined 
for the region Af < 0. 
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FIG. 32. The distribution of £ 2 values obtained from J/ipK + Monte Carlo samples equivalent 
to the data. A fit to the standard \ 2 distribution is shown by the solid curve and yields 12.7 degrees 
of freedom. The vertical line indicates the £ 2 value obtained from the J/ipK + data sample, and 
Monte Carlo samples with larger values of ( 2 are shaded. 
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